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PREFACE 


This  thesis  serves  primarily  as  an  introduction  to  the  analysis 
and  design  of  linear  shift-invariant  discrete-time  systems  possessing  a 
two-time-scale  property.  As  done  for  continuous  singularly  perturbed 
(two-time-scale)  systems,  we  hope  to  exploit  the  relative  "spread"  between 
eigenvalues  to  obtain  reduced-order  simulation  and  design.  Then,  we  can 
propose  computationally  efficient  algorithms  to  generate  the  control  law 
to  be  implemented  on  the  high  order  model. 

In  Section  1.2  of  the  next  chapter  we  give  basic  definitions  and 
requirements  for  discrete-time  systems  to  possess  the  two-time  scale  pro¬ 
perty. 

In  Chapter  2,  we  propose  a  general  class  of  discrete-time  systems 
that  will  satisfy  the  two-time-scale  requirement  if  a  parameter  is 
sufficiently  small.  Then,  we  outline  how  a  general  discrete  system  may  be 
"fitted"  to  our  model.  By  slightly  modifying  the  derivations  in  [6],  we 
derive  sufficiency  conditions  on  our  parameter  such  that  our  two  limiting 
cases,  denoted  "H"  and  "V"  satisfy  the  two-time-scale  property.  This  will 
enable  a  "slow-fast"  block  diagonalization  to  be  possible.  Finally,  we 
conclude  this  chapter  with  sufficiency  conditions  for  "slow-fast"  block 
diagonalization  of  our  general  form  and  with  some  comments  on  the  dual 
nature  of  our  transformations. 

Chapter  3  uses  the  recursive  property  of  discrete  time  systems  to 
show  how  truncated  versions  of  the  N-delayed  sequence  of  x^(k+l)  may  be  used 
to  approximate  the  systems  fast  modes  as  delays.  Also,  by  assuming  that 
A^|  exists  both  forward  and  reverse  recursions  are  possible.  This  enables 


V 


us  to  "reflect  ahead"  terms  in  the  series  to  determine  the  homogeneous 

structure  at  some  desired  future  iteration.  This  idea  is  applied  to  both 

homogeneous  and  nonhomogeneous  systems. 

2 

In  Chapter  4,  O(^)  and  0(p,  )  approximations  of  control  laws  are 
designed  for  the  slow  modes  of  near  limiting  cases.  Then,  high-accuracy 
reduced-order  models  for  eigenvalue  placement  are  proposed  for  each  limit¬ 
ing  case.  This  is  then  extended  to  our  general  model  and  an  efficient 
algorithm  is  proposed  for  use  in  a  multiprocessing  control  environment. 
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1.  INTRODUCTION  -  BASIC  CONCEPTS 
1.1.  Problem  Statement 

Methods  for  approximate  control  of  large  scale  systems  have 
received  a  great  deal  of  attention  in  recent  works.  Of  these  methods, 
aggregation,  disaggregation,  and  singular  perturbations  seem  to  be  the 
most  well  known  [1  ] .  The  analysis  and  design  of  continuous  linear 
singularly  perturbed  (two-time-scale)  systems  has  been  well  documented 
[2,3,4].  The  multiple-time-scale  property  of  these  systems  has  been  used 
in  arriving  at  reduced  order  models  and  control  laws  for  these  higher  order 
"stiff"  models . 

Until  recently,  just  about  all  the  work  done  on  systems  possessing 
a  multiple-time-scale  property  has  been  done  on  continuous  systems.  The 
area  of  discrete  two-time-scale  systems  has  received  little  attention,  and 
with  the  increase  in  computer-based  adaptive  control  schemes  in  a  time¬ 
sharing  environment,  the  need  for  a  thorough  discrete  analysis  and  a 
simplified  design  procedure  for  systems  possessing  this  property  becomes 
essential . 

The  basic  model  assumed  for  the  system  dynamics  is  a  set  of  N 
linear  shift-invariant  difference  equations  of  the  form 

x  (k+1 )  =  Ax(k)  +  Bu  (k)  (1.1) 

where  x  is  the  vector  state  of  the  system  belonging  to  an  N-dimensional 
vector  space,  u  is  the  vector  input  belonging  to  an  M-dimensional  vector 
space,  and  the  matrices  are  assumed  to  have  constant  coefficients  and  be  of 
compatible  size  with  the  above  vectors.  The  variable  k  is  integer  valued. 
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There  are  three  important  sources  of  discrete-time  models  [9]. 

One  is  for  the  approximation,  usually  for  the  purpose  of  digital  simulation, 
of  the  system  of  continuous  differential  equations 

x(t)  *  Ax (t)  +  Bu (t)  (tQSt<t^)  (1.2) 

by  corresponding  difference  equations.  For  example,  if  a  first 
difference  is  used  to  approximated  a  derivative  then  the  following  approxi¬ 
mation  results. 

x(tt+A)  =  [I+AA]x(t.)  +  BAu(t  )  (1.3) 

where  A  is  the  time-discretation  interval.  For  a  discrete  model  obtained 
in  this  way,  (I+AA)  is  generally  nonsingular  since  A  must  be  small  for  a 
good  approximation  to  a  continuous  time  system. 

The  application  of  Euler's  approximation  to  continuous  singularly 
perturbed  systems  has  been  investigated  in  [ 10]  and  has  produced 
interesting  results. 

Another  source  of  the  discrete-time  model  is  sampled-data  systems. 
In  sampled-data  systems,  a  continuous  system,  such  as  represented  by  (1.2), 
is  driven  by  an  input  specified  at  discrete-time  points  and  has  output  and 
state  variable  analytically  available  only  at  discrete-time  points. 

The  standard  example  of  a  sampled-data  system  follows  when  u(t) 
is  a  piecewise  constant  function  of  time,  i.e. 

u(t)  -  u(tt)  t1it<ti+1 

and  the  state  and  output  are  sampled  at  discrete  time  points  t. .  The 
sampled-data  discrete-time  model  can  then  be  computed  from  the  continuous 
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differential  equations  by  using  u(t)  =  uCt^),  then 


x(t)  =®(t)ti)xi+J  0(t,T)BdT  u(ti)  (ti^.t<tf).  (1.4) 

Ci 

By  letting  t  =  t^+^  we  obtain  our  discrete-time  model  with 

Ci+1 

A  =  *(ti+l*V*  B  =  /  nt1+1,T)BdT. 

t. 

1 

Since  the  transition  matrix  $(t,T)  is  nonsingular,  the  matrix  A  is  non¬ 
singular.  This  class  of  discrete  models  is  ideal  for  computer  control 
analysis  of  continuous  systems. 

The  final  class  of  discrete-time  systems  results  when  the  physical 
system  is  inherently  discrete.  Such  models  are  common  in  economic, 
biological,  and  sociological  systems.  The  system  matrix  A  for  this  type 
of  system  may  well  be  singular.  So,  the  problems  associated  with  large 
scale  systems  are  for  the  most  part,  preserved  for  these  discrete-time 
models . 


1.2.  Basic  Definitions 

In  this  section  we  define  what  is  meant  by  a  two-time-scale 
linear  shift-invariant  discrete  systems.  These  defining  properties  were 
originally  proposed  in  [8]  . 

Consider  the  linear  shift-invariant  discrete-time  system 


x1(k+l) 

x2(k+l) 

'll 


'21 


+ 

X' 

x 2  (k) 

- 1 

CM 

PQ 

_ 1 

u(k) 


(1.6) 


I 

I 


where 


i. 


;‘£>’pN  “<k)€R" 

ii.  A^j ,  Bj  i, j  *  1,2  are  real  constant  coefficient  matrices 
of  suitable  size 

iii.  The  eigenvalues  are  concentrated  in  two  distinct  groups 

of  N  and  N,  such  that 
s  f 

a.  N  *>  N  +  N, 
s  f 


b.  Denote  X  c  4  max  X  .  i  €  N, 

t  t  i  r 

X  4  min| X  I  j  6N 
s  j  j  s 


Xf<<Xs<<0(l/Xf) 

or,  the  maximum  possible  magnitude  of  any  eigenvalue 

of  our  system  must  be  within  some  neighborhood  of  1. 

Note :  By  defining  an  upper  bound  on  the  magnitude  of  our  eigenvalues  of 

0(  /X,)  ,  our  fast  modes  will  always  be  stable  and  only  the  slow  modes  can  be 

unstable.  This  property  will  be  assumed  throughout  the  text. 

Now,  given  (1.6)  with  conditions  i,  ii,  and  iii,  then  there  will 
N 

exist  a  basis  in  R  such  that  (1.6)  takes  the  form 


xg(k+l)  Ag 

-  —  »  —  —  *  _ 
xf(k+l)  0 


0  x  (k) 


1  A 
I  Af 


Xf(k)  |Bf 


+  ---  u(k) 


(1-7) 


where  i .  x  (k)  €  R 

3  Nf 

ii.  xf(k)  €R  . 

This  transformation  is  not  necessarily  modal.  However,  multiple  and  complex 


conjugate  eigenvalues  are  naturally  grouped  together  in  either  A  or  A,. 

s  f 
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I 

I 

I 

I 

I 

I 

I 
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Under  this  transformation  (1.6)  is  said  to  possess  the  two-time-scale 
property  if 


Then 


tnax|\  j  (Af )  | 


minjX.  (A  )|  . 

J  X.  O 


\f  «xs  «  0(l/\f). 

In  other  words,  if  there  is  sufficient  "gap"  between  the  eigenvalues 
A  and  Af,  then  (1.6)  possesses  the  two-time-scale  property. 

S  L 

We  can  express  this  desired  property  in  the  form  of  matrix 
norms,  since  for  any  nonsingular  matrix  S, 


of 


max|\i(S)|  <  ||s|| 

[min|\ j  (S )  1 1  _1  <;  'Is'1!! 
where  we  will  define  our  norm  as 


!|a||  =  [Xmav[A*A]]^  *  =  conjugate 

ulcLX 

transpose . 

Therefore , 

min|\.(A  )|  >  ||a  ^||  ^  (lower  bound) 
£  1  s  s 

max| Xj  (Af)|  £  ||Af  11  (upper  bound) 


and  our  requirement  for  a  two-time-scale  property  in  terms  of  these  bounds 
becomes 


Next,  we  will  investigate  a  general  syscem  structure  for  which 
there  exists  a  basis  such  that  decomposition  into  our  block  diagonal  form 
is  assured  providing  a  certain  parameter  is  sufficiently  small. 
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2.  A  CLASS  OF  TWO -TIME -SCALE  DISCRETE  SYSTEMS 
2.1.  Introduction 

In  this  chapter,  we  introduce  a  model  for  a  class  of  two-time¬ 
scale  discrete  systems.  We  then  give  a  simple  algorithm  that  will  attempt 
to  fit  any  discrete-time  system  to  our  general  model. 

In  the  remaining  sections  of  this  chapter  we  attempt  to  isolate 
the  slow  and  fast  subsystems  through  similarity  transformations.  The 
original  proof  of  existence  of  these  transformations  was  used  in  the  study 
of  continuous  singularly  perturbed  system.  However,  we  will  show  that  by 
modifying  the  derivation,  we  can  adapt  these  ideas  to  a  class  of  discrete 
systems.  Also,  we  can  define  two  unique  transformations,  one  giving 
computational  simplicity  to  the  slow  and  one  to  the  fast  subsystems.  These 
transformations  will  become  very  important  in  our  design  chapter. 

We  conclude  this  chapter  with  a  brief  discussion  on  degenerate 
systems  and  some  comments  on  the  dual  nature  of  the  transformations. 


2.2.  System  Forms 

A  general  class  of  linear  shift-invariant  discrete-time  systems 
possessing  a  two-time-scale  property  may  be  expressed  as 


x1(k+l)' 

1  (1- j )  * 

Au  ;  ^  J\2 

x:(k) 

x2  (k+1) 

M.jA21  ;  M-A22 

*2  (k) 

+ 


u(k> 


(2.1) 


where 
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Q  11  f 

i.  x^  (k)  €  R  ,  x2  (k)  6  R  ,  u(k)€Rm. 

A  A 

ii.  A.  ,  B  ,  i,m  =  l,2  are  real  constant  coefficient  matrices 
£m  m’ 

of  appropriate  order, 

iii.  0  <  j  <  1. 

iv.  p.  is  a  small  positive  parameter  such  that  0  <  p  <  1. 

Note:  Throughout  the  text,  the  following  notation  will  be  consistent 

hi  -"■*V 

For  j  =  1  and  j  = 0  we  define  two  limiting  cases  which  we  will  show  to  have 


dual  characteristics  to  one  another.  These  limiting  structures  are 


A 

Structure  "H"  *  ( j  =1 ) 

x^k+1)  Au  J  A12 

x2(k+l)  p.A21  J  PA^ 


xx(k) 
x2  (k) 


(2.2) 


Structure  "V"  ^  (>0) 

xi  (k+i)  A1X  ;  ha12 

x2(k+l)  A21  J  pA22 


(2.3) 


For  0  <  j  <  1  we  will  prove  that  the  general  form  satisfies  the 

two-time-scale  property  under  similar  sufficiency  conditions  as  either  of 

2 

the  limiting  cases.  In  this  analysis,  however,  when  0(p)  and  0(p  )  reduced 
order  models  are  considered,  only  limiting  or  near  limiting  cases  will  be  of 
practical  importance. 

In  the  next  section,  the  problem  of  putting  a  general  two-time¬ 
scale  discrete  system  into  our  general  form  is  considered. 
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2.3.  Permutation  of  States  and  Establishing  of  Parameter 

In  this  section,  we  outline  a  straightforward  procedure  for 
transforming  an  Nth  order  discrete-time  system  into  our  proposed  general 
form.  This  procedure  could  be  easily  programmed  for  easy  analysis  of  any 
discrete  system. 

Given  the  linear  shift  invariant  discrete-time  system 

x  (k+1)  “  Ax(k)  +  Bu(k)  (2.4) 

where  i.  x(k)€RN,  u(k)€RM 

ii.  A, B  are  real  constant  coefficient  matrices, 
iii.  There  exist  two  distinct  groups  of  open  loop  eigenvalues 
numbering  Ng  and  N^  such  that 

a.  N  +N,-N 

s  f 

b.  If  =  max  |  \  |  6  N, 

A  j  J 

\s  -  min|\i|i^Ns 
then  \f  «\s«  0(l/\f). 

Permute  the  states  finding  the  ratio  of  the  norm  of  the  homogeneous  part 

of  Ne  of  the  states  to  the  norm  of  the  homogeneous  part  of  the  remaining  N 
f  s 

states.  Continue  this  until  all  state  combinations  have  been  tried.  Retain 
the  permuted  order  such  that 

Ml  A 

*  minimum 
and  the  permuted  system  looks  like 


y(k+i) 

a* 

A  !  B 

*y(k) 

_2(k+l)_ 

C  J  D 

_z(k) 

(2.5) 
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where 


Ns  Nf 

y  (k)  6  R  ,  z(k)€R  . 


Once  this  minimum  has  been  £ound,  we  have  identified  our  p.  parameter 


as 


IIdII 

*  '  liar- 

By  evaluating  ,|B  ||  and  ||c  || ,  we  can  then  establish  whether  or  not 
our  general  model  can  fit  the  system.  We  now  discuss  some  different  cases 
when  evaluating  these  norms. 

i.  If  || Bi|,  ||c||  are  both  on  the  order  of  ||a||  or  greater,  then  the 
system  cannot  be  put  into  our  general  form. 

ii.  For  the  remaining  cases,  a  closer  examination  must  be  made.  First, 
define  the  constants 


Then  find 


(1-j)  ^  0 


j  =  9 


log  (Y) 
logGO 

log(p.) 


>  0 

>  0. 


We  now  can  give  a  general  system  form  which  has  our  class  of 
systems  as  a  special  case 


y(k+l)" 

•  9+1 

A  |  \i  B* 

"yOO 

z  (k+1) 

j  M-D* 

z(k) 

where 

mV  = 

B, 

p3C*  -  C, 

*—• 

fl 

* 

Q 

Now,  by  observing  the  sum  S+0 ,  we  can  cover  all  possible  cases 
Case  a:  1-e  £  8-*6  £  1-K  0  <  £  «  1 


or,  for  within  some  small  neighborhood  of  1. 


If  0-H3  satisfies  these  oounds ,  the  system  is  referred  to  as 
"normal,”  and  has  the  following  properties: 

1)  The  system  fits  our  model  with  9  ■  (1-j)  and  0  ^  (j). 

2)  The  new  subsystem  matrices  B*,  C*,  D*  have  all  been  normalized 
with  respect  to  A  and  therefore, 

IIaII  as  ||b*||  as  ||c*||  «  ||d*!|  . 

Case  b:  0  £  040  <  1-e 

With  these  bounds,  we  are  approaching  condition  i.  Therefore,  the  system 
usually  cannot  fit  the  model. 

Case  c:  1+e  <  0-f0  <  30 

In  this  case,  the  system  is  referred  to  as  "degenerate,"  and  has  the 
following  properties: 

1)  The  system  can  be  made  to  fit  our  model  by  appropriately  scaling 
down  x  and/or  y  such  that 

$*  =  0  -  60 

0*  =  9  -  &9 

until  0*+0*  satisfies  case  a,  then,  the  norm  of  B*  will  be  scaled 
down  by  p,^  and  C*  will  be  scaled  down  by  resulting  in 

||a||  ||d*||  >  1Ib*||,  ||c*||. 

2)  The  system  has  ''weaker"  coupling  between  subsystems  than  a 
"normal"  system  with  as  similar  value  of  p. 

As  shown  in  a  later  section,  a  degenerate  system  will  have  some 


advantages  over  its  corresponding  normal  form. 
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12 


| All  -  1-46,  UbII  -  1.11,  Hell  -  .173,  ||dI1  -  .1414 


|a| 


.09685 


at 


A 


.11849 


.7603 


«  4  a  -  l£2igd 
log[nJ 


(1-j)  -  9 


iog[y] 

login] 


.91362 

.11738 


9+S  ■  1.031  ss  1  (normal). 


Therefore  the  system  fits  our  model  and  in  this  form  looks  like 

*1(k+l) 

*2(k+l) 

x3(k+l) 

,x4  (k+1 ) 

But  notice,  ^  (•  H^38)  is  close  to  B.  Therefore,  to  approximate  this 
as  a  limiting  case,  ideal  for  perturbations  analysis,  let 


1.0 

-.2 


.6 

1.3 


(.11738) 


.5261 
,  0 


1.052 
.  921/ 


,(•91362) 


.844' 


-.844  .844, 


.5162 

0 


1.0325^ 

1.0325> 


xx(k) 
x2  (k) 
x3(k) 
Lx4(k) 


.11738 


ss 


1, 


.91362 


which  identifies  limiting  case  H  and  looks  like 


- 1 

X 

»-* 

'-'t 

1  .6  .4  .8 

x3(k) 

x2(k+l) 

-.2  1.3  0  .7 

x2  (k) 

x3(k+l) 

/  0  1.032\  / .5162 ,  1.0352\ 

x3  (k) 

f  M  ) 

_x4(k+l) 

\-1.032  1.032 j  \  0  ,  1.0352/ 

_x4(k)_ 

p.  -  .09685. 
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In  the  next  few  sections  the  value  of  p,  becomes  critical  in 
determining  whether  or  not  a  system  can  be  exactly  decomposed  into  slow- 
fast  subsystems.  If  conditions  are  not  met,  a  scaling  of  states  as  well 
as  a  permutation  may  be  necessary. 


2.4.  Slow-Fast  Block  Diagonalization  Using  Riccati  Type  Equations 

For  continuous  systems  possessing  a  two-time-scale  property,  [6] 
has  developed  sufficient  conditions  to  block  diagonalize  the  system  so  that 
one  subsystem  contains  all  Ng  slow  modes  and  the  other  subsystem  contains 
all  fast  modes.  The  sufficient  condition  is  based  on  the  norms  of  the 
original  subsystem  matrices.  By  slightly  modifying  the  derivation,  we  can 
apply  Lemma  1  of  [6]  to  achieve  sufficiency  conditions  needed  to  block 
diagonalize  both  structure  "H"  and  "V"  such  that  they  take  the  form 


xs(k+l) 

_As 

0 

xs(k)' 

+ 

B 

s 

xf (k+1) 

0 

Af 

xf  (k) 

,Bf. 

u(k) 


(2.6) 


and  the  submatrices  Ag  and  A^  satisfy  the  two-time-scale  property. 

Our  purpose  for  handling  only  these  two  limiting  structures  is 
so  that  we  can  first  understand  the  origin  of  our  dual  transformations. 

These  two  system  structures  motivate  straightforward  methods  of  solving  dual 
algebraic  Riccati  equations  needed  to  define  our  slow-fast  transformations. 
These  solutions  are  found  using  a  successive  approximations  scheme  where  it 
is  critical  what  we  select  as  an  initial  iterate.  These  ideas  will  carry 
over  to  the  next  section  where  these  transformations  become  special  cases  of 
our  general  class  (0£j<Ll). 
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Part  1 :  Block  diagonalization  of  form  "H" 
Given  r- 


(k+1) 

An  ; 

A12 

xx(k) 

x2 (k+1) 

.^21  ! 

P‘A22_ 

x2  (k) 

let 

y1(k)  =  x:(k) +Lx2(k)  (2.8) 

yi(k+l)  =  x^k+1)  +I^2(k+1) 

*  AllXl(k)  +A12x2(k)  +^LA2lX;L(k)  +  p,LA22x2(k) 

-  (A^LA^y^k)  +  (-AnL+A12^a22-^LA21L)x2(k). 

We  want  to  find  L  such  that 

-Al:L  +  A12  +  M.LA22  -^LA21L  =  0.  (2.9) 

In  [11]  a  unique  solution  to  this  equation  is  found  through  partitioned 
submatrices  of  the  modal  matrix.  If  only  an  approximation  of  L,  namely  Lq, 
is  available,  an  efficient  method  of  improving  this  initial  guess  can  be 
obtained  using  a  successive  approximation  scheme  on  the  algebraic  Riccati 
equation  (2.9)  with  Lq  as  the  first  iterate. 

In  our  case,  the  calculation  of  a  modal  matrix  is  something  we 
computationally  want  to  avoid.  In  general,  there  still  remains  the 
problem  of  picking  a  good  first  guess  Lq.  For  this  limiting  structure, 
however,  the  choice  is  obvious.  For  |4  small, 

-AUL  +  A12  +  Oflj.)  =  0 
Lo  3  A11A12’ 

If  A^  exists,  which  we  will  show  to  be  a  necessary  condition  for  (2.1)  to 
possess  a  two-time-scale  property,  we  can  define  the  following  successive 


approximation  scheme. 
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Letting  L  “  Lq+D  and  substituting  into  (2.9)  we  get 
-Au(Lo+D)  +A12  +  p.(Lo+D)A22  -M.(Lo+D)A2ia0-H))  =  0 

-A11D  +  hLqAo  +  M-DAq  -iaLqA21D  =  0  (2.10) 


where 


Ao  "  ^2  "  A21Lo 


-1 


and,  since  A^  exist,  we  naturally  define  the  successive  approximation 
scheme 

Dk+1  =  M,All[LoA21Dk  +  DkA21Dk‘LoAo‘DkAo]  ' 

If  we  define 

a=u.||Aotl,  b  =  |Ia21I|  IIlJI,  C  «  ||aJ[J|| 


M 


k  HlWJI 

we  can  apply  the  contraction  mapping  argument  of  [6]. 
Lemma  1 :  If 


c  <  j  (a-R.)'1 


or  if 


!u;}ll  <  i  0.IUJ!  ^I!a2.||  IIlJI)-1. 


‘21"  "-o' 

There  exists  a  unique  real  root  of  (2.10)  bounded  by 

-1 


0  <1  ||df||  £  2(a*>) 


0  £  !]  D||  <L 


2||AoH  ilLj 


HaJI  +Ha21||  IIlJI  ' 


(2.11) 


(2.12) 


Moreover,  this  root  is  an  asymptotically  stable  equilibrium  of  (2.11). 

Proof :  Appendix  A. 

Due  to  the  explicit  presence  of  n,  we  state  an  obvious  sufficient 
bound  on  p,  for  convergence  of  (2.11). 
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Corollary  2 :  To  put  structure  H  in  lower  block  triangular  form,  it  is 
sufficient  that  p,  satisfy  the  following  bounds 

0  £  >.  < - : - - - .  (2.13) 

3||a1-;!I(IIaJM|a21||  IIlJI) 

The  proof  follows  directly  from  equation  (2.12)  of  the  previous  lemma. 

To  complete  the  block  diagonalization,  let 
y2(k)  =*  x2(k)  -My1(k)  (2.14) 

y  2  (k+l )  =  x2  (k+1)  -  My1(k+1) 

°  ^21yl^  +  ^  (A^-A^L)^  00  -  M(A11  +p.LA21)y1(k) 

=  p.  (A22*A21L)y2  (k)  +  (p.  (A^-A^L^-MCA^-HJ-LA^)  +p,A21)y  1  (k)  . 

We  seek  M  such  that 

ti(^22'^21L)M"  M(A11  +^l^2i)  +^i  =  °-  (2.15) 

By  identifying 

As  =  11  (^22*A21L) 

As  =  An  +  p.ll21 

(2.15)  becomes 

A^  -  MAg  +m,A21  =  0.  (2.16) 

From  [14]  we  know  that  this  Lyapunov  equation  will  have  a  unique  solution 
since  Ag  and  A^  have  no  eigenvalues  in  common. 

Again,  through  successive  approximation,  we  can  show  that  M  is 
an  asymptotically  stable  equilibrium  of  the  difference  equation 

«k+l  ’  >*ti22Mk-i21U,k-Mk“21+i21lA[1l-  «-17> 


Lemma  2  of  [6]  again  use  a  contraction  mapping  argument  to  provide  a 
sufficient  condition  for  convergence  of  (2.17). 
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Lemma  3 :  If 


c  < 2  (a+b)  ^ 


or 


Kill  <i  dUJI  +IIXJ  IIlJI)-1 

then  (2.16)  has  a  unique  real  root  such  that  ]|m!|  =  0(p).  Moreover,  this 
root  is  an  asymptotically  stable  equilibrium  of  (2.17) 

Proof:  Appendix  A. 

Our  transformed  system  now  becomes 


71  (k+1 ) 

Ly2  (k+1). 


^All"t*J,LA21'* 


M-  CA22‘M'A21L')J 


Vj.00 

y2  00 

To  verify  the  two-time-scale  nature  of  this  system,  we  know 

Kll  =  HAll+iiLA21l| 

II Afl|  =  M-l|A22-A21Ll! . 

In  Chapter  1  we  showed  that  for  a  linear  shift-invariant  discrete-time 
system  to  possess  a  two-time-scale  property  it  was  sufficient  that 

KTl »  iia£!i  . 

In  our  case  we  need 

||(An+^LA21)'1|r1  »  ti||A22-&21Ll| 

which  will  be  true  for  every  p.  satisfying  Corollary  2.  Putting  this  trans¬ 
formation  into  matrix  form  we  have 


i 

/^\ 

rM 

V-/ 

r— 1 
>> 

• 

\  ;  l 

\(k) 

—  ”  “  ~ 

-  -  —  — 

y  2  (k> 

-M  J  I2-ML 

x2  (k) 

Ns  Nf 

where  y^(k)€R  ,  y2  (k)  €  R  ,  and  1^  and  I^  are  Ng  and  order  identity 


matrices  respectively.  This  nonsingular  transformation  has  the  attractive 
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feature  that  its  inverse  is  simply  given  as 


*1 (k)  I1-LM  |  -L  yL(k) 


x2  (k> 


M  i  I2  y2(k) 


Applying  this  transformation  to  our  nonhomogeneous  system  "H1 


y  ^  (k+i )  (a11+^la21) 

y2(k+l)J  [  0 


y1(k) 


Ha22-A21L)|  |_y2(k)  L-MB1+(I-KL)B2 


B1+LB2 


yL(0)  =  x:(0)  +  Lx2(0)  (2.19) 

y2(0)  =  -Mx1(0)  +  (I-ML)x2(0). 

Note  that  the  y^(k+l)  subsystem  is  totally  independent  of  the  M  trans¬ 
formation  matrix. 

Part  2  :  Block  diagonalization  of  form  "V" 


Given 


(k+1) 


A11  .  ^A12  Xl(k) 


x2(k+l)j  [A21  ,  [x2(k) 


y2(k)  =  x2(k)  +  Lx1(k) 


(2.20) 


y2 (k+1 )  =  x2 (k+1 )  +  Lx^ (k+1 ) 

=  A21x1(k)  +u,A22x2(k)  +  LA1]_x1(k)  +  p.LA12x2  (k) 

=  (p.A22+M-LA12)y2(k)  +  (A^-p-A^L  +  LAj_ ^p-LA^L^  (k) 


We  want  L  such  that 


A21  ~^22L  +  ^11  "m'LA12L  =  °‘ 


For  p.  small  a  good  first  approximation  would  be 
A21  +  LA^  +  0(p)  =0 

Lo  =  -A2 1A11* 


(2.21) 
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We  can  again  apply  the  contraction  mapping  argument  of  [6].  Letting 

L  =  L  +D 
o 

A21  "^A22(Lo+D)  +  a0+D)All '^(Lo+D)A12(Lo+D)  =  0 

-|iAoLo  -p.AQD  +  DAU  -^DA12Lo  -^DA12D  =  0  (2.23) 

where 

Ao  =  ^22  +  LoA12  * 

Since  A^  exists,  we  again  naturally  apply  the  successive  approximation 
scheme 

Dk+1  "  *  (AoLo  +AoDk  +DkA12Lo  +DkA12Dk)AU* 

If  we  define 


We  can  again  apply  Lemma  1  of  [6]. 

Lemma  4 :  If 

c  <  j  (a+b)"1 

mil  <  5“  (IIaJI  +||a12I1  IllJI)"1.  (2.24) 

Then  there  exists  a  unique  real  root  of  (2.22)  bounded  by 


such  that  D  is  an  asymptotically  stable  equilibrium  of  (2.23). 
Proof :  Appendix  A. 


We  now  give  the  obvious  bounds  on  p. . 
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Corollary  5:  To  transform  form  "V"  into  upper  block  diagonal  form,  it  is 
sufficient  that  p.  be  within  the  bounds  defined  by 

OSK  - - - - - -  .  (2.25) 

3iu;;n  (ikjMixji  iiLjb 

The  proof  follows  directly  from  (2.24)  of  Lemma  5.  To  complete  the  block 
diagonalization  let 

yx(k)  =  xx(k)  -  My2(k) 

y  L  Os+l )  =  xx  (k+l )  -  My 2  (k+1 ) 

=  [A11-p.A21L]y1(k)  +[  (A11-p.A21L)M-MCiA22-hj.LA12)  +uA12]y2(k). 

We  can  identify 

As  =  A11  '^l1, 

Af  =  ^(A22  +LA12)* 

Therefore,  we  want  M  such  that 

AsM  -  MAf  +  p.A12  =0.  (2.26) 

Again  from  [14]  we  know  that  this  Lyapunov  equation  will  have  a  unique 
solution  since  Ag  and  A^  have  no  eigenvalues  in  common.  Through  successive 
approximation,  we  can  show  that  M  is  an  asymptotically  stable  equilibrium  of 

“k+1  ■  *AU<il2'"Mk+Mki22  +VVi12)-  <2'27> 

Lemma  6 :  If 

c  <  j  (art)'1 

llAuli<i(l|Aj|H.||A12ll  IlLjl)-1. 

Then  the  root  to  (2.26)  is  an  asymptotically  stable  equilibrium  of  (2.27). 


L 
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Proof :  Appendix  A. 

Our  transformed  system  now  looks  like 


y:(R+i) 

Air^i2L 

!  o 

i 

^00 

y2(k+i) 

0 

i  M>  (^22+LA12^ 

y2  00 

and  again  it  is  obvious  that  our  two  time  scale  property  is  satisfied 
since 


(AU-UA^)'1!!’1  >5>  “11*22  +LX 


12" ' 


This  transformation  in  matrix  form  becomes 


yjds) 

IX-ML  ! 

-M 

xx(k) 

72  (k). 

L  ' 

x2(k) 

(2.28) 


N 


where  y1(k)€R  ,  y2(k)€R  and  1^,  I2  are  and  order  identity 

matrices  respectively.  This  also  has  a  simple  inverse  given  by 


x1(k) 

I,  !  M 

1^00 

_x2  ^ 

_-L  |I2-LM 

72 

Applying  this  transformation  to  our  original  nonhomogeneous  system  we  get 

u(k)  (2.29) 


y^k+i) 

s 

(Aii-p.Ai2L)  1  0 

1  l 

i  *  ' 

+ 

(I-ML)B1-MB2 

y2(k+D_ 

0  1  ^(A22+lA12)_ 

72  (k)_ 

LB1+B2 

where 

y1(0)  =  (I-ML)x1(0)  -  Mx, (0) 
y2(0)  =  Lx1(0)  +  x2(0). 

Note,  in  this  case  the  y2(k+l)  subsystem  is  independent  of  the  M  trans¬ 
formation  matrix.  This  is  dual  to  the  results  obtained  in  part  one,  and 
will  become  very  important  later  on. 
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2.5.  General  Case  -  Block  Transformation 

Up  to  now,  we  have  dealt  exclusively  with  our  limiting  case 
models  "H"  and  "V."  This  was  done  to  introduce  the  reader  to  the  two 
types  of  transformations  and  to  easily  see  the  formulation  of  the  matrix 
Riccati  recursion  equations  as  an  application  of  [6]  and  an  analog  to  what 
was  done  in  [5] .  What  we  do  now  is  to  apply  both  transformations  to  the 
general  model  and  not  to  isolate  limiting  cases.  This  will  enable  us  to 
develop  sufficiency  conditions  for  either  transformation  to  be  applicable 
to  the  general  model.  This  idea  will  become  extremely  important  when 
reduced  order  control  laws  are  considered. 

Since  many  of  the  derivations  are  similar  to  section  (2.4),  some  of 
the  repititious  steps  will  be  eliminated.  Given 


x^k+1)  _  Au  H(1'j)A12  *1(k) 

x2(k+l)  |^A2i  M-A22  x2  (k) 

where 

i.  All  matricies  have  constant  coefficients 

Ns  Nf 

ii.  x^(k)€R  ,  x2(k)€R 

iii.  0  <  j  £  1 


(2.30) 


iv.  p,  is  some  small  positive  parameter  <  1. 


Representing  p.  A^  by  A^ ,  let  us  consider  the  general  transformation 


y2(k)  =  x2(k)  +  Px^(k)  where  P  *>  L  as  j  =>  1. 


Our  notation  has  been  changed  due  to  our  convention  in  the  limiting  cases 
that  L  =  0(1)  and  M  =  0(p,). 

Even  in  this  general  case,  P  assumes  the  role  of  the-  independent 
transformation  matrix.  And,  as  a  result,  the  homogeneous  submatrix  and 
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control  submatrix  of  the  y2  00  (fast)  state  variables  will  be  dependent 
only  on  the  P  transformation  matrix.  We  will  denote  this  as  the  "F" 
transformation.  Continuing  the  derivation 

y2(k+l)  =*  x2(k+l)  +  Pxx  (k+1)  (2.31) 

=  [A21  +  PA11-A22P-PA12P]x1(k)  +[A22  +PA12]y2(k) 

and  obviously  we  want 

A21  +  PA11  '  A22P  '  PA12P  "  °  (2-32) 
which  is  equivalent  to  (2.21)  except  for  P  =  L.  Again,  we  would  like  to 
apply  the  contraction  mapping  argument  of  [6]  and  obtain  a  convergent 
successive  approximation  scheme  so  as  to  be  able  to  iterate  a  solution  to 
(2.32).  However,  an  initial  guess  to  such  an  iteration  is  not  as  obvious 
as  before. 

Make  the  substitutions 

A21  =  ^21’  A22  =  ^^22 *  Ai2  *  ^(1"j)^i2  (2-33) 

in  (2.32),  and  we  obtain 

^jA21  +  PAU  -^A22P  -  Ri(1‘j)A12P  =■  0. 


For  (j.  small  and  (0£j<l),  a  good  0(p.)  initial  guess  is  not  so  obvious. 
Let  us  try  a  general  form  of  Lq  used  in  (2.21) 

Po  =  ■p,jA21All’ 

Substituting  this  into  (2.32)  we  get 


',jA 
d+j) 


•Is 


p.  A21  +  (-p.  A21A11)A11 +p.A22p.  A2lAu-P  A21A11M.  J  A12A2 1^ ^ A11 
1 


A22A2 1A11 


A2 1A11A12A21A11 


-V; 

-i 


(2.34) 
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or  P  =  Pq  +  0(^(1+j)) 

*  -^Ja21Au  +  °^1+j>- 

Thus,  P-PQ  “  0(ji^  ^ )  and  Pq  will  always  be  an  order  of  ^  greater  than  any 

error  between  P  and  P  . 

o 

Now,  since  exists,  we  can  naturally  define  the  successive 
approximation  scheme  to  improve  our  solution  to  P  as 


Pk+1  =  (A22Pk  +PkA12Pk'A21)All 


where  PQ  is  specified  to  be  Thus,  if  we  can  show  the 

DaP-PQ  is  an  asymptotically  stable  equilibrium  or 

“kfl  ■  Vo+AoPk+DkA12Po+DkA12Dk>AU 


(2.35) 


(2.36) 


where 


Ao  =  A22  +  PoA12' 


We  can  guarantee  the  P  is  an  asymptotically  stable  equilibrium  of  (2.35) 
If  we  define 


a-||A0l!,  b  -  ||P0||  IU12II,  c-  dk 


‘Aoil  ilPoii 


We  can  again  apply  the  contraction  mapping  argument  of  [6].  Therefore,  the 


iteration  converges  when 


'Allil  <3  (11^2  +  PoA12ll  KJIj' 


c  <  ~  (a+b)-1 


(2.37) 


Making  substitutions  (2.33) 


lAt'Jll  <  J  (]]^2-k)A21A-1V1-jA12||+VJi21A[11||  II)*1. 


However,  for  every  matrix  H,  and  constant  u 
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M  «  \v\  ||h!I 

!|a"J||  <  j  (m.||a22-a21a^a12!!  h-m.!|A21a‘J;||  IIaJI)'1 


or,  M>  is  bounded  by 


0  1  M.  < 


jIIa^JII  (1|a22-a21a‘}'a 


111121 


+  ,l|A2iAir 


v12' 


Also,  D  will  be  bounded  by 


o  £  ||d| 


2/j)IIa0II  IIpJI 

<IIAoIHIPoII  11*1211 


and  our  upper  triangular  system  looks  like 


Xx  (k+l) 

3 

(A11~A12P)  !  A12 

xx(k) 

1 2  (k+1)_ 

.  o  ;  <aj2+pa12)_ 

y2  GO 

(2.38) 


(2.39) 


(2.40) 


To  complete  the  block  diagonalization 

yL(k)  -  Xl(k)  -Qy2(k)  (2.41) 

y L (k+l )  =  (A11-A12P)y1(k)  +  [  (A11-A12P)Q-Q(A22 +PA12) +A21]y2(k). 

We  seek  Q  such  that 

(AirA12P)Q'Q(A22+PA12)  +A21  =  °*  (2.42) 

We  can  identify 

As  =  *11 'Hi* 

Af  "  *22  +  P*12  ’ 

Thus,  we  seek  Q  such  that 

AsQ  -QAf  +  A^  -  0.  (2.43) 

Again,  this  Lyapunov  equation  will  have  a  unique  solution  since  Ag  and  A^ 
have  no  eigenvalues  in  common.  However,  we  can  show  that  Q  will  be  an 
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asymptotically  stable  equilibrium  of  the  matrix  difference  equation 

Qk+1  ■^5<A12P<>k  +  ‘’k<422+PA12>-A12)'  <2'44> 

Using  Lemma  2  of  [6]  we  know  that  this  equation  will  converge  to  a  stable 
equilibrium  if 


c  <  \  (a+b)'1 


which  implies 


<liA22+Poi12:l  +!'A12:I  »Po’I> 


-1 


which  is  automatically  satisfied  by  (2.38).  Making  the  substitution 


Aij  Aij  in  (2‘44) 


*k+l 


4  ^  ^All(A12Qk  +  Qk^A22  +  PA12^"A12') ' 


Therefore 


Q  ^  OOj.1'^)  -  M 


P  ^  OftJ)  ^  L 


j-0 

j  *0 ' 


(2.45) 


If  we  had  used  the  independent  block  transformation  on  the  x^(k)  subsystem 
first,  (call  this  the  "S"  transformation) 


yx(k)  =  xx(k)  +  Px2(k) 
y2(k)  «  x2(k)  -  Qy ^  (k ) . 


(2.46) 


A  similar  set  of  results  will  follow.  We  will  outline  these  results  to 
avoid  too  many  repetitious  derivations.  The  two  matrix  recursion  equations 
are  given  as 


Pk+1  *  A11(A12 +PkA22_PkA21Pk') 

Po  “  A11A12 

i1,  Qk+1  "  (A22Qk'A21PQk‘QkPA21 +A21)A11 
„  _  .  -1 


(2.47) 


(2.48) 
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Both  of  these  nonlinear  matrix  difference  equations  converge  to  stable 
equilibrium  values  if 

IIa'JII  <  j  (IIa22-a21pJ| +||p0||  IIa^II)"1. 

Showing  p.  parameters  explicitly, 

IIaJJI  <  j  (||p.A22-p.jA21A^H(1‘-j)A12H  +||lJ.jA21l|  ,,|A‘Jp.<'1*j)A12|!)'1 

ilA-ll  <  (||a22-a21a^a12!I  +||a21||  ||a^a121I)'1 

which  gives  us  the  following  bound  on  p. 

0  £  \i.  <  - r - r - - - t -  .  (2.49) 

||A^Ji|(||A22-A21A^jA12il  +!|a21II  ||a[Ja12||) 

Also,  it  is  easily  seen  that 

P  =  0(jj,1_j)  Q  =  0(p.j).  (2.50) 

This  similarity  in  bounds  for  transformations  "F"  and  "S"  is  very  important. 
Both  of  these  transformations  have  been  derived  using  our  general  two-time¬ 
scale  form.  Therefore,  if  p.  satisfies  both  bounds,  then  either  transformation 
can  be  used.  To  avoid  confusion,  in  future  sections  the  superscript  s or  f 
will  be  used  in  conjunction  with  the  P  and  Q  transformation  matrices  to 
distinguish  which  transformation  was  used. 

2.6.  Degenerate  Cases 

We  defined  a  "degenerate"  case  as  one  where 

1-K  <  8+p  <  »  0  <  e  «  1.  (2.51) 


Let  us  take  the  degenerate  model 
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y(k+l)‘ 

A  |  +®B 

y(k) 

z(k+l) 

"0"  r  "  - 
1-*  C  i  m»d 

z(k) 

(2.52) 


where  9  and  p  are  bounded  by  (2.51). 

9  9 

Let  us  implicitly  carry  p,  ,  p.“,  and  p,  with  3,  C,  and  D  respectively 
Then,  the  sufficient  condition  for  block  diagonalization  (transformation  F) 
becomes 


With  explicit  p.  parameters  we  have 

IIa_11|  <  j  (Hm.d-u.0ca“ V9 b1!  +||^0ca‘1||  lip-8 b'i|  )~ 1 
IIa"1||  <  ^  (||d-p. (9+6“1)ca"1b||  +p.9+?'1||ca'1|!  bib"1. 

Therefore 

0  *  3|(a-1||(|I^'9+»-1^-1b||+,^-1||ca-1!I  || i{|)  ' 

However,  if  the  system  where  "normal,"  the  bound  on  p.  would  be 

3|1a-1||  (||d-ca-1b||  +IICA"1!!  bll) 

We  now  state  a  lemma  enabling  us  to  show  that  the  upper  bound  on  p.  for  a 
degenerate  case  will  always  be  greater  than  for  a  "normal"  case.  This  will 
enable  us  to  concentrate  on  "normal  cases,"  realizing  that  by  modeling  the 
degenerate  system  as  "normal"  we  will  always  achieve  results  in  approximation 
convergence,  etc.,  that  are  equal  to  or  better  than  the  "normal”  case. 

Lemma  7 :  If  a  discrete-time  system  can  fit  our  model  only  in  a  degenerate 
form,  such  that 

l‘«  9+0  <  “ 


(2.53) 


(2.54) 
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then  the  upper  bound  on  p,  needed  for  bj.ock  diagonalization  is  greater  than 
the  corresponding  normal  case.  In  other  words, 


3||A'1||(lb->*(8',9'l)0A'18||  ^<S-'«-1)||cA-1||  ||bII) 
^  1 


3|1a-1!1  (||d-cab||  +i!cA"i||  ;  |  b|  I  ) 


111  il_!l 


(2.55) 


Proof :  (2.55)  implies 

!Id-cab||  4-  Hca-1!!  ||b||  >  ||D-vi.e+0"1CA-1Bj|  +  p.04S'1||ca'11|  ||b|| 
||d-ca_1b||  -  ||d-u(94^'1)ca'1b!|  >  Cx9+e"1-i)||cA"1!|  ||bI|. 

Let 

f  =  ||d-ca~lb|J  -||d-u(9‘^“1)ca'1b|| 


and  let  us  examine  the  bounds  of  F  as  e+p  "swings"  through  the  degenerate 
range.  Then  our  inequalities  become 

as  9+9  -  1  F  -  0  >  O04?'1-l)i|cA‘1!|  ||bI|  (2.58) 

as  9+9  -  »  F  -  |[ca_1b||  >  (P.9‘^-1-1)|'|CA"1||  ||b(|.  (2.59) 

Thus  F  is  bounded  such  that 

o  <  f  <  ||ca"1b|{ 


as  long  as  1  <  9+9  <  <=,  which  defines  a  degenerate  system.  Therefore, 
(2.55)  is  to  be  satisfied,  then  both  upper  and  lower  bounds  on  F  must 
satisfy  (2.56).  Looking  at  (2.58) 

0  >  Oi(e+S'1)-l)|(cA~1||  Ml 

0  >  H  (e4€_1)  -  1 
i 


if 


(2.60) 
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which  will  be  true  for  every  degenerate  value  of  9+S  since  —  0<^«  1. 

Just  for  completeness  (2.59)  gives 

||ca"1b||  <  (p,946"1-i)I|ca‘1||  ||bI| 

;1ca'L°II  +  1 

II CA’1!!  ||b|| 

which  is  always  true  if  (2.60)  is  true.  Therefore  (2.55)  holds. 

The  impr ctance  of  this  is  that  we  may  consider  our  general  model 
as  an  "upper  bound"  for  many  degenerate  systems  which  also  satisfy  a  two- 
time-scale  property.  Also,  since  many  of  the  models  and  algorithms 
proposed  in  this  thesis  are  (J.  dependent  in  terms  of  accuracy  and  convergence 
properties,  higher  values  of  u  can  be  used  with  degenerate  systems  and  still 
obtain  good  results.  However,  even  more  important,  by  modeling  a  degenerate 
system  in  our  general  form  (i.e.  B  and  C  are  no  longer  normalized  to  A),  we 
can  expect  better  results  (convergence  of  algorithms,  reduced  order 
approximations,  etc.)  than  for  a  similar  system  (same  |i  parameters)  with 
"normal"  submatrices  B  and  C. 

We  can  conclude  this  section  by  stating  that  by  performing  our 
analysis  and  design  procedures  on  our  "normal"  general  form,  we  can  expect 
any  degenerate  form  of  this  normal  form  to  have  more  accurate  results. 

2.7.  Duality  of  Transformations  "S"  and  "F" 

In  two  preceding  sections  sufficiency  conditions  for  the  block 


diagonalization  of  our  general  two-time-scale  model  and  for  the  two  limiting 
cases  were  carried  out  separately. 


The  purpose  of  this  was  to  point  out  the  dual  nature  of  the 


transformations  that  were  derived.  In  particular,  transformation  S  may  be 
considered  the  "slow"  transformation  since  it  will  allow  slow  mode  modifi¬ 
cation  independent  of  the  iteration  (as  shown  in  a  later  chapter) . 
Likewise,  F  may  be  considered  the  "fast"  transformation  for  the  similar 
reason.  This  property  will  become  very  important  when  we  consider  eigen¬ 
value  P  lac  erne  nt . 

When  attempting  to  block  diagonalize  the  general  model,  for  either 
transformation  we  found  very  similar  bounds  on  p  for  every  j  including  our 
limiting  cases.  What  this  means  is  that  for  p  small  enough,  either  trans¬ 
formation  may  be  used.  Thus,  we  can  choose  whether  we  want  our  slow  or 
fast  state  vector  to  be  independent  of  Q^. 

Our  two  transformations  in  general  form  look  like 


where 


‘*l(k)] 

=  S 

\  GO 

x-  (k) 

y,GO 

1—  Z. 

—  Z.  - 

r  _ 

1 

x1(k) 

y:(k) 

x2  (k) 

*  p: 

1 2 (k). 

Qf  -  0<]u(1'j)) 


ri-Qsps : 


S  = 


P 

L.  S 


F  = 


-Qf  !  I-QfPf 


Pf  =  0<pj) 


(2.64) 


(2.65) 


QS  =  O0iJ) 


Ps  =  0(p.<:i‘j)) 


and  the  transformations  are  related  by 


S  -  <F-V. 


In  block  diagonalizing  our  limiting  cases,  P1  -  L-O(l),  i=s,f  and 

Q.  =  M^0(p).  This  was  done  to  preserve  the  order  of  L  and  M.  However,  for 
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0 <  j  <  1 ,  the  order  of  the  matrices  as  well  as  the  off  diagonal  system 
matrices  will  vary.  When  these  transformations  are  independently  applied 
to  our  general  two-time-scale  model,  the  resulting  nonhomogeneous  slow-fast 
systems  are : 
using  "F" 

xa(k+l)  (Ai;L-)iA12Pf)  !  0  xg(k) 

xf(k+l)_  _  0  1  ^(A22+LA12)  |  xf  (k) 

(I-Qfpf)B  -QfB, 

+ - f  *  -  -  -  u(k) 

p  b1+b2 

using  "S" 


f  s 

where  x^(k)  is  independent  of  Q  and  (k)  is  independent  of  Q  . 


3.  RECURSIVE  PROPERTY  AND  DELAY  APPROXIMATION 


3.1.  Introdution 

This  chapter  serves  primarily  as  an  introduction  to  reduced  order 
modeling  of  discrete  two- time-scale  systems  fitting  our  general  model.  By 
taking  advantage  of  the  recursive  property  of  discrete- time  systems,  we 
will  be  able  to  show  how  the  x^Ck+l)  and  (k+1)  subsystems  may  be  expressed 
as  Nth  order  matrix  difference  equations  in  increasing  powers  of  H. 

x1(k+l)  =M1xi  +  p,M2xi(k-l)  +p.NMNxi(k-N) 

(k+1)  =K]xi  +  |J.K2xi(k-l)  +^\Nxi(k-N) 

By  truncating  this  "expansion"  after,  for  example,  two  delays 

2 

(0(|i  )),  we  will  establish  our  first  attempt  at  a  reduced  order  "slow" 
model . 

Then,  since  it  is  assumed  that  A^  exists,  we  can  take  advantage 
of  the  backward  recursion  for  the  x^(k+l)  subsystem, 

xL(k)  =A^x^(k+l)  "Aj^p.^  ^A12x2(k)  (3.1) 

to  "reflect"  ahead  the  model  of  our  system  at  some  future  iteration.  This 
will  result  in  an  intuitively  satisfying  derivation  of  our  slow  subsystem 
found  through  transformation  in  Chapter  2.  By  continued  application  of  the 
backwards  recursion  of  x^(k)  to  our  Nth  order  difference  equations,  the 
existence  of  our  "P"  transformation  matrix  will  become  apparent  as  we  con¬ 
verge  to  a  "steady  state"  or  "slow"  model. 

We  conclude  this  chapter  with  a  look  at  the  nonhomo geneous  case 
to  see  how  the  control  variables  fit  into  this  derivation. 
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Before  going  on  to  the  next  section,  let  us  define  a  term  that 
will  be  used  throughout  this  chapter. 

Definition 

A  matrix  difference  equation  such  that  the  next  state,  x^(k+l)  is 
dependent  on  the  previous  N  time  instances  such  that 

xi(k+l)  = f (Xj (k) ,  Xj(k-l),...,  Xj(k-N))  (3.2) 

where  j  is  not  necessarily  equal  to  i,  will  be  refered  to  as  an  N-delay 
sequence  of  x^(k+l)  in  x^(k-j),  j  =  0, 1,2  , .  .  .N. 

3.2.  Expansion  of  x^(k+l)  Subsystem  into  an  N-Delay  Sequence 

A  linear  shift-invariant  discrete-time  system  is  by  nature  a  set 
of  N  1-delay  sequences.  Thus,  using  the  recursive  property  of  discrete 
time  systems,  it  can  be  expressed  as  a  function  of  the  last  N  events,  a  N- 
delay  sequence.  In  this  straightforward  set  of  substitutions,  we  hope  to 
show  how,  as  more  delay  terms  are  added,  there  contribution  to  the  next 
state  x^(k+l)  decreases  as  a  power  series  in  H-. 

In  the  following,  we  will  assume  that  the  system  is  in  our 
general  linear  shift-invariant  discrete  two-time-scale  form.  Also,  for 
ease  of  notation,  we  will  use  the  A^j  submatrices  until  the  substitution 
A^j  is  needed. 

We  now  hope  to  show  that  for  any  system  that  can  be  put  into  our 

N 

general  two-time  scale  form,  the  x^(k)  states  can  be  approximated  to  0(p.  ) 
by  an  autonomous  (function  of  only  delayed  versions  of  itself)  delay-N 


sequence . 
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Given 


where 


x1(tefl) 


A12  xl^k> 


(3.3) 


x2(k+l)  A2l  !  A2?  x2(k) 


N  Nf 

i.  x2(k)  £R  S  ,  x2(k)  €  R 

ii.  Aj_j’  3“^*^  are  rea^  constant  coefficient  matrices 
satisfying  our  general  two -time- scale  structure 


x^k+1)  =A11x1(k)  +A12x2(k) 

—  A^  ^ x^  (k)  4*  A^2^2 (k- 1 )  4"  A12A22 ^2  ^k~  ^  ^ 


continued  substitution  of 


x2(k-j)  =A2lx1(k-(j+l))  +^22x2  (k-l(j-t-l)) 


(3.4) 


results  in 


x^(k+l)  =  A11x^(k)  +a^2A21X1  +  A12A22A21Xl^k~2^+* . .  .+ 

A12A22 ^2 1X1 (k"N^  +  A12A22^X2 ^k‘N  ^ 


XA 

Now,  making  the  substitutions  A„  =  p.  A^  , 

x1(k+l)  *A11x1(k)  +  p,<'L'j)A2lx1(k-l) 

+  n(i'1)A12^A22(iJA2lx1(k-2) 

+  ...+  p,<'i‘j)A12(|iA22)N_  VA21xL(k-N) 


+  ^,(i-j)^12(^A22)  x2(k_N) 


(3.5) 
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“  A11X1  (k>  +  ^i2^2lXl^k'L^  +  M‘2A12^22^2lXl^k'2^ 

N  ^  A  W  —  1  A  M-i- 1  —  1  A  A  M 

+  . . .+  ^  A12A22iN  1A21x1  (k-N)  +  t*W+i  JA12A22i  x2  (k-N)  (3.6) 

The  significance  of  (3.6)  is  that  for  small  perturbations  in  x2 (0) ,  the 

entire  dynamics  of  the  x^(k)  subsystem  may  be  represented  accurately  to 
N 

00*  )  by  an  N-delay  sequence  in  x^(k). 

Thus,  by  using  the  process  of  state  augmentation,  we  can  simulate 

the  x^(k)  states  to  a  desired  accuracy. 

2 

For  example,  if  we  want  0(p  )  accuracy,  our  model  would  look  like 
XL (k+1)  *AL1xL(k)  +pA12A21x1(k-l) 


Let  2(k)  *  x^(k-l) 


xL(k+l) 

z(k+l) 


A11 

^12A21 

x^(k) 

I 

0 

^  1 

v-' 

eg 

_ *J 

(3.7) 


where  x^(0)=x^q  ,  z(0)*0. 

For  Ng«Np,  this  model  is  good  for  p,  sufficiently  small.  However, 
for  or  for  p  too  la*se,  the  order  of  the  model  can  surpass  the  actual 

system  model  to  achieve  reasonable  accuracy.  As  a  first  attempt,  this  model 
proved  to  be  very  important  as  will  be  seen  in  future  sections. 


3.3.  Use  of  Backward  Recursion  to  Establish  Steady  State  Models 

In  the  previous  section,  the  backward  recursion  (3.1)  of  x^(k) 
was  never  used  in  obtaining  our  Nth  order  difference  equation  for  x^(k+l) 
(3.6).  In  this  section,  however,  the  general  form  of  (3.1), 


xL(k-j)  “A’Jxj^k-j+l)  -A^^1  j:>A12x2(k-j) 

-J-0,1,2, .  (3.8) 

is  applied  once,  and  then  N  times  at  each  iteration  to  obtain  "s 
state"  models  of  the  system  at  the  next  or  Nth  iteration. 

Again,  given  system  (3.3), 

x1(k+l)  =A11x1(k) +A12x2  (k) 

x2(k)  =A21x1(k-l)  +A22x2(k-1) 

(3.10)  into  (3.9)  gives 

x^k+1)  =A11x1(k)  +A12A2lx1(k-l)  +  Ai2A22X2  ^k_2') 
using  (3.8)  for  j  = 1,  into  (3.11) 

x1(k+l)  =A11x1(k)  +A12A2lA1^x1(k) 

+  A^2  (A22”A2lAllA12^X2 

—  (A^  +  A^2A2^A^)x^  (k) 

■V  A^2  (A-22  *^2  2  )  ^2  ^  ^ 

x2(k-l)  =A21x1(k-2)  +A92x2(k-2) 

*  [A22‘A21A11A12^X2  (k"2)  +  A2lAUXl (k"  L) 
xL(k+l)  =  (A11+A12A21AL^)x1(k) 


now,  into  (3.12) 
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+  (A12(A21-A21A11A12)A2lA1L)x]L(k-l) 

+  (A12(A22-A2lA^Al2)2)x2(k-2)  (3.13) 

Continue  these  substitutions  in  the  same  manner  using  the  general  equation 
x2  (k-j)  s  (A22  ™A2 ^A^ 2^A^2  ) x2  (k*  j  *  1 )  ■+•  A2  jAj^x^  (k~  j ) 

j  *0,1, .  (3.14) 


We  generate  the  following  delay-N  sequence. 


-1. 


xL(k+l)  =  (AI1+A12A21A11)x1(k) 


d"  (A12  (A22 -A2 1A11A12 )A2 2_A^^) x^(k-l)+.  • . . 
+  (A12(A22-A2lA^Al2)N"1A2lA^)x1(k(N-l)) 


+  (A12(A22-A21A1^A12)N)x2 (k-N) 


(3.15) 


Again,  making  the  substitutions  A _  =»  p,  A„ 


x^(k+l)  -  ( A L L  p. ( L  j)A12pjA21A1^)x1(k) 

+  ^(I'j)A12(^22"^J^2lAn^(1'j)^12)^21All)xl(k'1) 


+  ••••  + 


n-n*  *  i*  -1  (1-i)*  N-liA  -1 

+  (p."  j;A12(pA22-^JA21AupV  J;A12)1  pJA2lA1J)x(k-(N-l)) 

n-iV  *  i  a  -1  ci-  m 

+  (pU  ai2(^22'^  A2lAll>  j;A12)N)x2(k-N)  (3.16) 


since  p,  ^  ^  are  positive  scalers,  they  can  be  combined  to  give 


Xl(k+1)  =  (A11+nA12A21A“^)x1(k) 

+  M1  (A12  (A22-A21A11A12)A21A11)x1(k-l) 


N  *  A  A  -  1  A  M-  1  A 

+  ^  (A12(A22-A21A11A12)  A21A1J)x1(k-(N-l)) 

+  p.N+1_j  (A12  (A22-A21A^A12)N)x2  (k-N)  (3.17) 


In  arriving  at  (3.17)  it  is  important  to  understand  that  we  have 

applied  the  backward  recursion  (3.8)  once  for  each  delayed  term  in  (3.6). 

This  substitution  has  "reflected  ahead"  the  influence  of  the  x^(k-M)C^ 

th 

term  on  the  x^(k-(M-l))  term.  Conceptually,  you  may  look  at  (3.17)  as 
our  N-delay  sequence  model  of  our  two-time-scale  system  for  k>l. 

In  comparing  N-delay  sequence  models  (3.6)  and  (3.17),  we  notice 
that  every  term  on  the  right  hand  side  of  (3.6)  is  of  the  form 


Flf(xL-M)  M>  1 


where,  the  equivalent  delayed  term  of  (3.17)  is  given  as 


M+l 

k™g(Xl-M)  M>1 


The  homogeneous  part  of  (3.17)  (M  =  0)  has  0(1)  and  0(|A)  term  whereas  for 
(3.6),  the  homogeneous  part  possesses  only  an  0(1)  term. 


In  Chapter  2,  through  transformation  we  arrived  at  a  slow  sub¬ 
system  matrix  As .  For  the  "F"  transformation  this  matrix  was  composed  of 


As  "All  "  A12Pk 
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•k+1 


(A22Pk  +  PkA12Pk  '  A21)A11  Po  *  'A21A11 


For  k  =  0 


A 

s 


‘A11  +  A12A21A11 


and  using  A 


ij 


'AU  +  ^i2^21All 


(3.18) 


Notice  that  (3.18)  is  just  the  homogeneous  part  of  (3.17).  So, 
by  substituting  the  backwards  recursion  into  (3.6)  at  each  delayed  tern 
(M>1),  we  have  generated  the  model  of  the  system  at  the  next  iteration. 
Since  x^(k+l)  does  represent  our  slow  state  variable,  the  equivalence  of 
the  homogeneous  part  of  (3.17)  and  equation  (3.18)  seem  logical. 

It  is  of  interest  now  to  examine  truncated  models  of  (3.6)  and 
(3.17)  so  as  to  determine  which  would  be  more  a  accurate  model  of  the  x^(k) 
states.  Three  criteria  must  be  examined  in  making  a  judgement. 

1)  Number  of  matrix  multiplications  required  to  generate 
the  model 

2)  Order  of  model 

3)  Magnitude  of  4 


For  an  equal  number  of  multiplications,  (3.6)  looks  like 

x^k+1)  -A11x1(k)  +  MA^A^x^k-l)  (3.19) 


and  (3.17)  like 
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xL(k+l)  =  [Au  +  4A12A2lA1^]x1(k)  (3.20) 

However,  (3.19)  is  twice  the  order  of  (3.20).  For  small  H, 
(3.19)  and  (3.20)  will  give  approximately  the  same  accuracy.  Therefore, 
since  lower  order  is  usually  preferred  (3.20)  would  be  the  logical  choice. 

As  M1  becomes  larger,  the  0(P>)  terms  in  (3.19)  and  (3.20)  become 
more  significant.  Then,  the  way  (3.20)  is  derived  from  (3.19)  becomes 
important.  In  other  words. 

Given 

xx(k+l)  =AuxL(k)  +nA12A2lxL(k-l) 

use 

x1(k-l)  =A‘Jxx(k)  -n(L'j)A‘jA12x2(k-l) 

This  gives, 

x1(k+l)  =A11x1(k)  +pA12A2l[A1|x1(k)  -p/1  (k- 1)  ] 

=  [Au  +  4Ai2A2iA1^]Xi(k)  -  42  ^^A^A^^^  (k- l) 

(3.21) 

Equation  (3.20)  neglects  the  0(p.  ^)  term  which  could  be  signi¬ 
ficant  for  larger  p,.  When  this  is  the  case,  (3.19)  will  give  a  more 
accurate  model.  This  method  of  modeling  will  be  refered  to  as  "delay 
compensation",  since  we  are  using  delays  to  compensate  for  nonneglible 
fast  modes  (larger  p.) .  More  will  be  said  about  this. 

Now,  using  this  same  line  of  thought,  let  us  apply  (3.8)  N  times 
at  each  delay  term  in  (3.6)  for  M>1. 

If  the  backward  recursion  (3.8)  is  repeatedly  applied  to  (3.6) 
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N 

until  f  has  0 (p.  )  accuracy  where, 

x^k+1)  =  fQXj^Ck)  +  fjX^k-1)  +  . .  .+  fjjXj^k-M)  +  gNx2  (k-  (N+l) 

Then  f  =  A  +0(p.N) 

o  s 

and  fL  Af2  A  ....  AfN  =  0(pN) 

To  properly  show  this  using  induction  is  straightforward  but 
very  tedious. 

Here,  we  offer  a  simple  heuristic  explanation  mere’y  expanding 
the  series  to  N  *  4.  Given  our  original  system  equations, 


x1(k+l)  -AI1x1(k)  +A12x2(k) 
x2(k+l)  = A21x1(k) + A22x2 (k) 
using  the  following  substitutions, 

x2  (k-j)  =  (A22-A2lA^A12)x2(k-(j-l))  +  A2lAilxl(k_J) 
x1(k-j)  = Al[x1(k-(j+l))  -  a1|a12x2 (k-j)V. 

for  N  *  1, 

xL(k+l)  =  [A11+A12A2lA1^]x1(k)  +A12[A22-A21AL^A12]x2(k-l) 
»  [As  +  0(p)]x1(k)  +0(p)x2(k-l) 


for  N  =  4, 
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Xj/k+l)  “  fAll  +  A12A21All  +  A12(A22  *  A21A11A12)A21A11 
A12(A22  'A21A11A12^A21A11A12A21A11  +  A12('A22  'A2lAllA12) 

A21All^Xl(k) 

A12^A22  'A21A11A12^  A2lAUA12A21AllXl(k_1) 


-  Al2  ^A22  ”A21A11A12^  A21A11A12A21A11X1  <‘k'2) 


j  +  A12^A22  "  A21A11A12^A21A11A12A21A11A12A21A11  j 

3  1 

"  +  Ai2fA22  "A21A11A12^  A21AU 


x1(k-3) 


+  0(p.4) 


(3.22) 


Which  is  of  the  form 


x^k+1)  =  fgXj^k)  +  fIx]L(k-l)  +  f2x1(k-2)  +  f-jX^k-3) 


where 

fo  *A11  "  A12Pk  =As  +  ^ 

and 

Pk+1  "  (PkA12Pk  +  A22Pk  "  A21 

-1 

Po  =  ‘A21A11 

and  from  (3.22)  it  is  obvious  that  f^,  f^,  f^  and  g^A  0(p,4). 

Continuing  in  this  manner  results  in  the  P  transformation  matrix  recursion 

k 

of  Chapter  2.  Conditions  for  its  convergence  have  already  been  studied. 
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Remark 

The  recursive  substitutions  were  made  in  an  attempt  to  converge 
to  the  F  transformation  matrices.  By  performing  the  dual  substitutions, 
the  S  transformation  matrices  could  have  been  easily  obtained. 

Thus,  we  have  shown  that  by  applying  the  backward  recursion 
equation  (3.8)  repeatedly  to  the  N-delay  sequence  of  x^(k+l)  (3.6),  we 
approach  a  "steady-state"  or  slow  model  of  the  system  at  some  future  itera¬ 
tion  N.  We  identified  this  with  our  slow  subsystem  developed  through 
transformation  in  Chapter  2.  For  p,  small,  this  "steady-state"  or  slow 
model  will  prove  to  be  very  accurate  in  approximating  the  x-^(k)  states. 

For  p.  larger,  the  initial  deviation  of  x^(k)  due  to  fast  modes  may  be 
significant.  "Delay  compensation"  may  then  be  used  to  more  closely  model 
these  initial  deviations.  We  refer  this  initial  period  of  deviation  as  the 
"boundary  layer"[2]  .  The  existence  of  this  "layer"  will  become  apparent 
when  we  next  discuss  the  x^(k+l)  subsystem. 

3.4.  x2(k+l)  Expansion  and  the  Boundary  Layer 

Up  to  now,  there  has  been  no  attempt  to  simulate  the  x2(k+l) 
subsystem.  If  the  same  autonomous  expansion  is  attempted  such  that 

x2(k+l)  =  f  (x2(k) ,  x2(k-l),...  x^(k-N)) 


we  get. 
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N- 1  N 

4“ >a21A11  ^12x2  1  ^ i  i  (k~^) 


21111 


making  the  substitution  A, .  = a  A. . 
6  tj  J-J 


x2(k+l)  -  |xA22x2(k)  +  ^A21^^1  ^A12x2(k-1) 

*  I 

+  M.(-l'j')A2lA;LL)ijA12x2(k-2)  +  .  .  . 


i*  N-l  ( 1- i ) A 
d-^A^A^  j;A12x2(k-N) 


4-  ^A21A^1x1(k-N) 


(3.23) 


x2  (k+1)  =  ^22X2  ^  +  pA21A12X2  <‘k‘I') 


H^2iAllA,12X2  • 


a  N-1A 

+  pA2iAli  A12x2(k-N) 


4-  ^A^A^x^k-N) 


(3.24) 


Examining  (3.24),  we  see  the  exists  no  truncated  delay  sequence 
which  could  accurately  model  the  x2 (k)  states.  If  fact,  the  x^(k-N)  term 
will  always  be  dominant.  Thus,  the  lowest  order  simulation  model  for  the 
^(k)  states  is  just  the  original  state  equation 

x2  (k4-l)  =|j,jA21x1(k)  4-pA22x2(k) 

which  is  obviously  nonautonomous  and  therefore  does  us  no  good. 

Now,  if  we  try  to  derive  the  N-delay  sequence  for  x2(k4-l)  such 


^  o 


x0  (k+l)  -  f  (x^k)  ,  x^k-l), - -  x^(k-N)  ) 

we  ge  t 

x.,(k+l)  -A.^x^k)  +  A,.,A.,1x1  (k-L) 

2  N 

+  A,.,A  ,  iX^  (k-2 )+.  .  .+  A  ,.,A.,  ^x^  (k-N) 

+  A^+^x,  (k-N) 

Substituting  A,1-|i'iA,l>  A^-pA,.,,  A^,  -pl'L  ^A^ 

x2(k+l)  -  M*'*A2lx1(k)  +  M,L+'1A,.,A21xl(k-l) 

2+ 1 A  A  N+  i  '  * 

+  H.  JA00A21x1  (k-2)+. .  .+tx  JA,.,A,1x1  (k-N) 

N+tA 

+  |J.  A^x^k-N)  (3.25) 

Again,  wo  have  a  case  similar  Co  the  x^(k-t-l)  subsystem.  We 

N+j 

can  CoCaLly  simulute  Che  x.,  (k+L)  subsyscem  Co  0(p.  )  accuracy  using  a 

truncated  N-delay  sequence.  For  N>  l,  using  (3.25)  would  result  in  "delay 
compensation"  corms  implimenced  using  state  augmentation  which  has  been 
already  discussed. 

One  very  noticible  difference  between  the  x.,  (k+l>  sequence  and 
the  x^(k+l)  sequence  is  the  preservation  of  Initial  conditions.  When  wo 
truncate  (3.25)  we  eliminate  any  dependence  of  x.,(k+l)  on  \,(0).  Thus, 
from  k«0  to  k-k  ,  (3.25)  may  differ  considerable  from  t  no  ictuul  x,(k*l) 
states.  Since  all  fast  modes  are  stable  in  this  analysis,  the  question  is 
not  whether  the  truncated  version  of  (3.25),  x.,(k+l),  will  converge  to 
x  ^ ( k+ 1 )  ,  but  how  soon 
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I  *2  (k+1)  -  (k+L)  I  <  e  for  some  e>  0 

and  any  k>  0 

•k 

The  period  of  time  [0,k  ]  is  refered  to  as  the  "boundry  Layer"  in  the 
analysis  of  continuous  systems  possessing  a  similar  property  [2,4] .  For 

•k 

p  small,  k  can  be  as  small  as  1.  In  this  case,  only  a  few  terms  of 

(3.25)  would  be  needed.  However,  as  p  becomes  larger,  an  impractical 

number  of  terms  would  be  needed  for  simulation. 

Now,  let  us  apply  the  backwards  recursion  of  (3.6)  to  (3.25),  so 
N+i 

as  to  get  0(p  )  accuracy  for  each  delay  term.  After  considerable 

algebraic  manipulation,  (3.25)  becomes  of  the  form 


where 


x2  (k+1)  =  f0xL(k)  +  f  ixi L)  +  f2xi^k'2)+-  •  •+fNx1(k-N) 


fo  ~  (A21A11  +A22A21A11  "  A2lAllA12A21All  +  *  *  *  ^ 


(AU+A12A2iAi1+  .  .  .  • ) 


(3.26) 


"■(PN)(AS)5Cl(k) 


fLAf2^f3--^fN  =  0^N+j> 


where  PxT  is  the  Nth  iteration  of 
N 


Pk+1  ~  (PkAl2Pk  +  A22Pk  ‘  A21->AL1 


-1 


Po“'A2lAll  k  =  0>L> - N 


and 
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AS  =A11  '  A12Pk 

The  similarity  between  these  expansions  and  the  "F"  transformation 
should  now  be  completely  clear.  This  transformation  is  given  again  for 
convenience 

(k)  =  xg  (k)  +  QFxf  (k) 
x2(k)  = -pFxs(k)  +  (I  -  pFQF)xf(k) 

By  applying  the  backwards  recursion  (3.8)  on  (3.6)  and  (3.25),  we 

are  finding  our  system  model  at  some  future  iteration  N.  In  a  two-time-scale 

system,  the  farther  N  is  away  from  k  =  0,  the  smaller  the  x^(k)  state  vector 

N 

will  be.  In  fact,  it  will  have  decayed  to  0(p,  )  proportion. 

This  therefore  verifies, 

xx(k+l)  =Agxs(k)  +  0(tiN) 

x9(k+l)  =  -PFAsxs(k)  +0(nN)  k>N 

for  our  N-delay  sequences. 

So,  for  p.  sufficiently  small,  we  can  approximate  x^(k+l)  by 

xL(k+l)  =AsxL(k)  (3.27) 

and  the  x2(k+l)  states  as  a  function  of  x^(k+l)  states  ?k>  0 

x2(k+l)  =  -PFAsxs(k)  =  -P^k+l)  (3.28) 

or  just  an  output  variable. 


However,  if  we  use  (3.27)  and  (3.28)  and  p.  is  large,  these 
approximations  may  deviate  considerably  during  the  boundry  layer.  To 
improve  the  boundary  layer  approximations ,  delay  compensation  terms  may  be 
added  (i.e.  x^(k-j)).  In  this  case,  the  delay  term  with  maximum  j  terminates 
all  backwards  recursion  substitutions  instead  of  the  x^(k)  term. 

In  the  next  section  we  will  give  an  example.  Then,  in  the  follow¬ 
ing  section  we  will  see  how  the  control  variable  u(k)  would  enter  into  our 
N-delay  sequence  for  x^(k+l)  states.  We  will  limit  this  discussion  to  the 
x^(k+l)  subsystem  since  the  "steady-state"  dynamic  of  X£(k+1)  states  will  be 
determined  by  the  x-^(k-f-l)  states. 


3.5.  Example  of  Reduced  Order  Modeling 


The  model  of  a  steam  power  system  discussed  in  [19]  is  given  as 


-2  0  0  0  -4 

4.75  -5  0  0  o 

z  =  0  .167  -.167  002 

0  0  2  -2  0 

0  .025  .02333  .035  -.1125 

To  put  this  system  into  singularly  perturbed  form  we  first  re¬ 
index  the  state  variables  using  the  permutation 


P1  =  ^e4’  e5  *  e2  ’  e3  ’  el-) 

Where  e^  is  an  elementary  column  vector  whose  ith  entry  is  1. 


With  y  =  Pz,  the  new  system  becomes 
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•.1125 

0 

0 

-4 

0 


.02333 

.167 

2 

0 

0 


.035 

0 

-2 

0 

0 


0 

0 

0 

-2 

4.75 


.025 

.167 

0 

0 

-5 


Now,  we  scale  the  states  using  x  =  Sy,  where  S  =  diag(4,l,  1,2, 1) .  The  re¬ 
sulting  system  is 

l_Tll25  .0932  .14  0  .1 

0  -.167  0  0  .167 

0  2-200 
-2  0  0  -2  0 

000  2.375  -5 

which  is  of  the  proper  form  in  accordance  with  [5] . 

If  this  system  is  sampled  at  .8  second  intervals,  the  discrete 

model  is 


x(k+l) 


.9014 

.  1179 

.0525 

.0167 

.02104 

-.0196 

.8743 

0 

.025 

.02934 

-.0071 

.7342 

.20175 

.013 

.021067 

-.75 

-.0557 

-.032 

.19357 

-.014076 

-.306 

-.01694 

-.011 

.14278 

.013217 

x(k) 


with  eigenvalues 
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Xx  2  -  . 8777  +  . 1054 j 

X3  =  .0179 

X4  5  =  .2055  ± ,0236j 

therefore  N  =2.  N ,  =  3 
s  *  f 

I1auII  =  .9415,  !Ia12 li  =  .0625,  A21  =  .8184 
ilA22il  =  .2441 
The  system  is  degenerate  with 


|j,  =  .259253 


which  is  less  than  its  defined  upper  bound  of  .27158.  Using  one  iteration 

,  f  .  . 

of  our  P  successive  approximation 


P 


f 

1 


-.025827 

1.0207 


-1.029 

-.09667 


.4817  -.06101 


we  arrive  at  our  slow  subsystem 


xs (k+1) 


.875521 

-.059322 


.1748599 


.8779818  j 


x„(k) 

a 


with  eigenvalues 

kl  *87675  +  .  10184 j 


J 
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We  now  wish  to  model  the  entire  system  dynamics  using. 

^(k+1)  =Asx1(k)  xL(k)SR2 

Where  the  remaining  three  states  appear  as  output  functions  or 

x2  (k)  =  -P^  (k)  x2  (k)  €  R3 

To  compare  this  second  order  model  with  our  fifth  order  model,  the  response 
of  the  system  to  an  initial  perturbation  of 

xQ-(l,  -.8,  .5,  .2,  .6) 

is  plotted  versus  its  lower  order  approximation  (Figures  3.1  to  3.5).  The 
actual  (high  order)  states  will  be  designated  x^(k),  while  the  approximated 
state  will  be  identified  by  x^^Ck). 

„  f  ,  , 

As  you  can  see,  using  just  one  iteration  of  result  in  very 
good  accuracy. 

3.6.  Nonhomo geneous  Case:  Delay-N  Sequence 
Given : 


x^k+1)  *A11x1(k)  +A12x2  (k)  +  BjU(k) 
x2(k+l)  =A21xL(k)  +  a22x2  +  B2u  (k) 


(3.29) 


where  our  parameters  are  again  assumed  implicit. 

In  this  section  we  will  simplify  the  derivations  of  section  3.1 
and  briefly  describe  how  the  appearance  of  a  control  u(k)  will  fit  in  to 
the  delay-N  sequence  expansion  of  a  general  two-time-scale  discrete  system. 
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Iterations,  (k) 


FP-S379 


Figure  3.4.  Reduced  order  approximation  of  state  x,  (k)  . 


From  3.29 
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*2  (k)  *A2lxl(k'l>  +A22x2(k-1)  +B2u(k-1) 

x^(k- 1)  =A1^x1(k)  -  AL^Al2x2(k-l)  -A^B^Ck-l) 

x2  (k)  — A^^A^^x^(k)  ~  A2  ]_A^^A^2x2  (k~  1)  -  A2^A^^B^u(k-l) 

+  A22X2  +  B2u(k-1) 

*  CA22‘A2lAUA12^X2  <‘k_1)  +  [B2'A2lAllBl^u<;k_1) 

“t”  ^2 1^11X1^^ 

x2  (k+L)  =A1]x1(k)  +Al2[A22-A21A‘^A12Jx2(k-l) 

+  A12[  (B2-A91A”^B1)]u(k-l)  +A12A21A1^xl(k)  +  BjU(k) 

Xx (k+1)  =  [A11+A12A21A1^]x1(k)  +  Al2  [A22'A21A11A12^X2  (k'1'1 
+  Al2[B2-A91A1^BL]u(k-1)  +  B^u  (k) 

By  successive  substitution  of 

x2 (k-j)  =  [A22-A21A1^A12]x2(k-(j+l)) +  [B2-A21A1^B1]u(k-(j+l)> 

+  A2lAUXl(k-:i)  (3.30) 

we  again  are  "reflectintng  back"  or  shifting  H  terms  one  unit  forward  and 
thus  approximating  a  "future"  system  with  less  system  order.  As  a  result 


we  get 
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x^k+1)  =  (A11+Al2A2lA1^]x1(k)  +  B^k) 

+  Ai2  tA22’A'21AllA'12^A2 lAllX]/k”^ 

+  A12^B2'A2lALlBl^U(k-1) 

-1  2  -1 

*+•  A-^2  ^22*^2 1^11^12  ^  ^2 1^1 1X1  ^“2 ) 

+  A12^A22_A21A11A12^  1fB2~A21AUBl^u(:k"2) 
+ . + 


oil 


N  *  a  *  _  ].  •»  N-l'  -1 

+  ^  (a->2"A''  lAllAP  ]  Ao iA  1 ^  “  (N “  ^  ^ 

M-l*  *  a  _ia  .(N-2)  *  *  -l 

+  ^  Al2^A22*A2lALlAl'1  IB->  ”A2lA  i_i.B|_]uC*c-  (N- i)  ^ 


N+L-  t A  *  A  - 1  •*  M 
+  M-  JA12[A^1-A  ,1AuA1  ,]  Xa  fk-N ’l 


N*  »  *  - 1'  i  N  -  O  '  *  -1 

+  V-  A 1 2  [ A2 2 _A2 1A L  LA 1  ^ 1  [B2"A21AUBllu(k*N) 


With  this  expansion,  we  can  now  detine  some  low  order  models  for  the  \  ^  (.  k  +  L ) 
subsystem.  To  0 (p.") 

To  Ofp.) 


x  ^  lk+ 1)  “  A  j  ^x  ^  (JO  +  B  (.ki 


(.3.  331 


To  0(n  ) 


-1, 


xL(k+l>  -  [ALL  +  (j-A1aA.,1A1  ^x^k)  +  B^u(k) 


+  ^12  [®2”A2  iAi_j|BlluCk- L> 


(.3.34) 


Ln  matrix  form 


x  L(.k+l) 

zfk-t-l) 

(A  ^ 12A  ’lALl'  '  2  ^  1  ”A  ’ lA i  i®  i ^ 

. j . -  - . 

o  :  o 


xL(k) 


2(k) 


l 


;  u  t  k ) 


L  i 

_  J 


f3.35) 

This  model  possesses  some  interesting  properties  some  of  which  will  be  nude 


apparent  in  the  remaining  chapters. 


4.  REDUCED  ORDER  CONTROL  PROBLEMS 


4.1.  Introduction 

In  the  previous  two  chapters  we  have  shown  the  existence  of  the 
slow-fast  phenomena  for  any  system  fitting  our  general  model.  Up  to  now, 
our  reduced  order  modeling  has  been  applied  only  to  system  simulation.  In 
this  chapter  we  begin  the  problem  of  designing  control  laws  based  on  our 
reduced  order  subsystems  and  applying  these  control  laws  to  our  higher  order 
mode  Is . 

When  designing  stabilizing  control  laws  for  linear  shift-invariant 
discrete -time  systems,  perhaps  the  most  frequently  used  methods  are  the 
optimal  linear  regulator  design  and  pole  placement.  This  chapter  will  deal 
with  reduced  order  pole  placement  procedures.  Reduced  order  regulator 
designs  is  another  research  topic  in  itself  and  will  not  be  considered  here. 

When  we  consider  any  reduced  or  approximated  procedure,  the  two 
important  factors  are : 

1)  computations  involved 

2)  accuracy  when  implemented  on  high  order  model. 

The  previous  chapters  helped  to  show  us  the  computations  involved 
in  attaining  some  reduced  order  models.  Since  we  have  assumed  stable  fast 
modes,  we  took  advantage  of  the  fast  "transient"  to  develop  "steady  state" 
models  of  our  original  system.  These  models  helped  to  verify  our  slow-fast 
block  diagonalization  transformations  and  to  introduce  the  "boundary  layer" 
phenomena  of  the  x9  (k)  subsystem.  In  this  chapter,  we  take  this  idea  one 
step  further  by  now  incorporating  a  feedback  control  variable  into  the 
design.  Since  we  will  be  dealing  with  "slow"  and  "fast"  eigenvalue  placement 
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we  will,  for  the  most  part,  find  it  easier  working  with  the  transformed  state 
vectors  xg (k)  and  x^(k)  rather  than  with  our  original  state  vectors  x^  (k) 


and  x2 (k) 


In  the  next  section  we  will  give  some  basic  differences  between 


singular  perturbations  in  continuous  systems  and  reduced  order  modeling  in 
two-time-scale  discrete  systems. 


4.2.  Correspondence  Between  Discrete  Two-Time-Scale  Systems  and 
Continuous  Singular  Perturbations 

Our  dual  slow-fast  transformations,  in  general  form,  are  given  as 


Xl(k)  =  (I-PSQS)xs0O  -  PSxf  (k) 
x2(k)  =  QSxg(k)  +  xf  (k) 


(4.1) 


xx  (k)  =  xg  (k)  +  Q  xf (k) 

x2(k)  =  -P  xg(k)  +  (I-PfQf)xf(k), 


(4.2) 


Each  of  our  original  subsystem  variables  is  a  function  of  our  transformed 
slow  and  fast  variables. 

In  many  instances,  when  a  design  model  is  being  formulated,  the 

fast  modes  are  neglected.  It  is  the  purpose  of  this  chapter  to  examine 

reduced  order  approximations  to  our  original  systems  when  the  fast  modes  are 

neglected  or  when  their  control  design  is  to  be  carried  out  separately. 

This  has  been  well  documented  for  the  continuous  case  [4,5]  where 

the  perturbation  made  in  the  original  system  is  singular,  and  the  order 

reduction  achieved  is  from  N— N  • 

s 
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However,  in  trying  to  adapt  these  ideas  to  the  discrete  case,  a 
fundamental  concept  between  two-time-scale  discrete  and  continuous  systems 
must  be  understood. 

In  continuous  systems,  the  perturbation  allowing  the  real  part  of 
the  fast  eigenvalues  to  go  to  -“  is  the  same  perturbation  involved  in 
stating  that  the  x^Ct)  subsystem  goes  to  equilibrium  instantaneously. 

Using  the  block  diagonalization  transformation  defined  in  [12]  on 
a  continuous  two-time-scale  system,  we  get 


X 

> 

O 

_ 1 

X 

s 

s 

0  1 

s 

>f 

_°  ;  h.2 

xf 

Neglecting  the  fast  modes  is  equivalent  to  letting  the 

ReUtA^]]  _*"co  or  xf(c)  =  0  Vt  >0. 


When  our  two-time-scale  system  is  put  in  singularly  perturbed  form,  it  looks 
like 


_*1 

A  1 
A12 

*1 

>  1 
s  • 

— 1 

A22_ 

.X2. 

By  letting  p  —  0  and  assuming  Re{\ [A22] ] < ^ ,  we  are  again  assuming  that 

Re[\[A22]}  -  “  ®  either  way,  the  perturbation  is  consistent.  In  the  discrete 

case,  the  fast  modes  are  not  totally  contained  in  the  submatrix  or  any 

other  submatrix  to  any  order.  Therefore,  letting  \[A^]  -0  and 

[X2  (k+l)-X2  (k)  ]  -0  are  not  equivalent. 

The  problem  is  that  for  the  general  structure  we  have  assumed,  the 
2 

fast  submatrix  to  0(p  )  is  given  by: 


A£ 


'22 


A21A11A12 


0(P2) 


and  it  is  therefore  very  difficult  to  isolate  the  fast  modes  in  the  original 
system.  This  is  possible  only  through  a  canonical  form. 

Another  more  serious  problem  of  discrete  two-time-scale  systems 
is  that  of  "relative  order"  of  system  parameters. 

What  this  means  is  that  the  p. -parameters  present  in  the  discrete 
model  also  correspond  to  the  order  of  the  fast  eigenvalues.  So,  for  small 
p, ,  an  independent,  reduced  order  "slow"  design  will  give  good  results, 
however,  since  our  fast  \'s  are  all  ready  0(p,),  an  accurate  reduced  order 
fast  design  is  almost  impossible  without  including  many  series  terms.  As 
a  result,  we  will  show  that  by  iterating  the  matrix  we  will  be  able  to 
solve  our  fast  problem  to  a  desired  accuracy. 

In  the  next  section  we  comment  on  controllability  aspects  of  our 
reduced  order  models. 


4.3.  Controllability  Aspects  of  Reduced  Order  Models 

In  this  section,  we  discuss  some  controllability  bounds  on  our 
two-time-scale  system  such  that  our  reduced  order  design  models  may  be  used. 
Given  a  system  in  our  general  form 


xx(k+l)]  _  |~  An  |  ^(1'j)A12]  [xx(k) 

--  -  f  “ 

*2(k+l)J  [y  A2 1  I  J  [x2(k) 


(A. 5) 


Satisfying  the  usual  requirements,  we  have  shown  that  for  p.  sufficiently 
small,  there  exists  a  nonsingular  transformation  T,  such  that 
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1' 
L~2  1 


-1 

xs(k) 

=  T 

xf  (k) 

where  =  0(1) 


xg(k+l) 

As0i)  !  0 

B>> 

xf(k+l) 

o  ;  Afoo 

_xf(k)_ 

+ 

Bf  (i-)_ 

u(k). 


(4.6) 


If  T  is,  for  example,  the  "F"  transformation. 


As  A11  "  ^  ^  A12PN 

Af  =  ^2  +  ^  (  j)]?NA12 

Bs  =  C!-Vn>B1“V2 
Bf  =  PNB1  +  B2 


where 


Pk+1  "  j)pkA12Pk +w'^22PknJ'jA21] A11 

Po  =  -4jA2  iAn 


(4.7) 


(4.8) 


Qk+1  *  All^  j)pNQk+,,lC!kA22 +^jQkPNA12"'i<'  ^A12^ 


(4.9) 


Q  = 

%  ^  A1 1a12 

As  the  P  and  Q  matrix  recursions  approach  equilibrium,  since  \(As0i))  and 
\  (A£ O-i*))  are  sufficiently  far  apart,  controllability  of  our  original  system 
pair  (A,B)  guarantees  controllability  of  our  subsystem  pairs 

(AsG0,Bs<O),  <AfGO,BfOO). 

However,  in  real  time  applications,  computation  time  requirements 
will  make  exact  equilibrium  values  of  P  and  Q  unreachable  in  general.  Thus, 

P  -  P„  +  0&iN+j) 

0  •  %  +  0<^N+1'b 
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resulting  in 

Ag  (n> )  =  As  +  0(nN) 

AfGO  =  Af  +  0(iiN) 

Bg  (l-t )  =  Bg  +  OfcN+1“J) 

BfCi)  -  Bf  +  04i N+j). 

When  this  is  the  case,  bounds  on  the  controllability  of  our  reduced  order 

models  are  best  defined  in  terms  of  "weak"  and  "strong"  controllability 

conditions  of  perturbed  systems  defined  by  Chow  [17].  For  a  given  finite 

number  of  iterations  N  of  the  P  and  Q  matrices,  it  is  necessary,  for 

separate  design  of  slow  and  fast  subsystems  that  our  original  system  be 

N 

strongly  controllable  in  p,  (i.e.  if  we  truncate  (4.8)  and  (4.9)  after  N 
iterations,  (4.5)  and  (4.6)  will  still  span  the  same  controllable  subspace). 
This  is  analogous  to  the  condition  defined  by  Chow  [ 18]  for  continuous 
singularly  perturbed  systems. 


Expressing  (4.6)  as 


xs(k+l)  _  As+ooO 
x  f  (k+1 )  0 


x  (k)  5  +  0(p,N+1_j) 


Af  +04OJ  [xf00j  [J3 f  +  0 4i  J ) 


We  can  readily  apply  Theorem  1  of  [12]  to  establish  bounds  on 
N 

p.  to  verify  the  strong  controllability  of  our  original  system. 


Theorem  7 : 

If  exists  and  if 

rank[Bg,AgBg, . . 

4N-la  , 

* ,As  “a1 

=  N 

s 

rank[ , . . 

■jN- 1*  , 

=  Nf 

then  there 

★ 

exists  p.  5(4. 10) 

is  controllable  V(i€(0,jj,  ] 

N 

In  other  words,  (4.5)  is  strongly  controllable  in  that  as  n  -0, 
our  fast  and  slow  subsystems  maintain  their  controllability. 

By  checking  the  conditions  of  Theorem  7,  we  can,  in  theory, 
establish  a  lower  bound  on  the  number  of  iterations  of  P  and  Q. 

2 

4.4.  0(h),  PC.  )  Control  Approximations 

The  remainder  of  this  chapter  deals  with  finding  control  laws 
based  on  our  reduced  order  models  and  observing  their  accuracy  when  imple¬ 
mented  on  the  higher  order  system. 

Control  law  design  of  continuous  singularly  perturbed  systems  has 
been  well  documented  [4,5,7].  The  reduced  order  control  law  designs  for 
eigenvalue  placement  [5]  and  the  optimal  linear  regulator  [4,7]  achieve 
0(p.)  accuracy  with  respect  to  exact  feedback  gains.  This  accuracy  is 
usually  very  acceptable  since  the  open  loop  fast  eigenvalues  have  real  parts 
that  are  0(1/h)  and  the  slow  eigenvalues  have  real  parts  that  are  0(1). 

Also,  when  expressed  in  block  diagonal  slow-fast  form  [12],  the  continuous 
singularly  perturbed  system  has  the  nice  property  that  !B^  ^  0(1/h)  while 

j  |  B  II  -  0(1).  This  enables  you,  to  0(h)  accuracy,  to  decompose  the  control 
into  independent  slow  and  fast  as  well  as  the  homogeneous  system  dynamics. 
Thus,  the  three  orders  of  magnitude  present  in  continuous  singularly 
perturbed  systems  (0(h),  0(1),  0(1/h)),  make  independent  0(h)  accurate 
approximations  and  decisions  very  acceptable  and  appealing,  to  large  scale 
systems  possessing  this  property. 


In  discrete  two-time-scale  systems  such  as  we  have  defined,  the 
"relative  order"  of  parameters  does  not  permit  such  0(h)  approximations  to 


be  acceptable.  The  | A^i!  =  OGi)  so  any  0(p.)  approximation  would  only  be 
meaningful  if  we  just  neglect  the  fast  modes  and  design  a  reduced  order 
slow  control  (|| A  1  =0(1)).  Likewise,  depending  on  p.,  0(p.^)  approximations 
could  be  unacceptable  for  both  fast  and  slow  subsystems.  Also,  in  the  model 
we  have  used,  no  assumptions  have  been  made  on  the  relative  order  between 
,  and  ■ ! B^J)  -  This  enables  us  to  encompass  a  greater  class  of  discrete¬ 
time  systems  which  have  open- loop  two-time-scale  properties.  Granted,  by 
making  order  assumptions  on  |b^!!  ,  independent  control  decomposition 

may  be  possible,  as  will  be  the  case  for  sampled-data  singularly  perturbed 
systems.  However,  for  now,  we  are  considering  purely  discrete-time  systems, 
so  the  proposed  generality  will  remain. 

As  a  result  of  these  problems  concerning  two-time-scale  discrete 

systems,  we  divide  our  reduced  order  control  law  design  into  two  parts. 

2 

In  this  section,  0 (p. )  and  0(p.  )  approximations  are  considered.  In  other 

words,  0  or  1  iteration  is  to  be  used  in  our  successive  approximation 

scheme  to  keep  computations  to  a  minimum.  Thus,  we  achieve  a  reduced  order 

2 

slow  representation  of  our  high  order  system  accurate  to  0(p.)  or  0(p,  )  for 
the  purpose  of  slow  control  law  design.  Of  course,  the  price  we  pay  is 
insufficient  accuracy  to  use  a  reduced  order  fast  model  to  generate  a  fast 
control  law.  The  only  requirement  we  can  pose  on  such  a  design,  is  to  place 
order  N  slow  mode  where  we  want,  and  keep  original  N-  fast  modes  as 

S  £ 

unaltered  as  possible. 

In  the  next  two  sections,  we  take  advantage  of  our  simple  successive 

approximation  matrix  recursions  to  obtain  fast  and  slow  reduced  order  sub- 
N 

systems  to  0(p  )  accuracy,  where  N  will  be  the  number  of  iterations. 
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This  enables  us  to  generate  highly  accurate  slow  and  fast  models  for  the 
purpose  of  separate  fast  and  slow  control  law  designs.  The  price  you  pay 
is,  of  course,  more  computations  than  the  0  4a)  approach,  however,  we  can 
now  specify  both  fast  and/or  slow  modes  to  a  desired  accuracy,  using  much 
less  computation  than  a  high-order  pole  placement  algorithm  applied  to  such 
an  ill-conditioned  system.  Also,  if  the  feedback  gains  are  to  be  altered 
in  some  adaptive  scheme,  this  reduced  order  method  offers  even  further 
computational  simplicity. 

2 

We  now  proceed  with  our  04a),  0(jj,  )  designs.  Before  considering 
040  approximations  to  our  general  model,  let  us  consider  one  of  our  limiting 
cases,  say  "V."  Here,  the  perturbation  p.  -  0  has  meaning  in  both  the  actual 
and  transformed  state  space.  Given  the  system 


x  ^  (k+1 ) 

_ 

x2(k+l) 

A11  |  M‘A1' 


x:(k) 

+ 

V 

x2(k) 

.y 

u(k) 


(4.11) 


N, 


where 


i.  x:(k)  €R  S,  x2(k)  €R  \  u(k)6Rm 


ii.  A^ j , Bj  are  real  constant-coefficient  matrices 

iii.  p.  satisfies  Corollary  5. 

Then  we  showed  that  there  exists  a  nonsingular  transformation  (2.28)  such 
that  (4.11)  may  be  transformed  into 


u(k)  (4.12) 


xg(k+l)‘ 

All"^A12L 

0 

xs 

+ 

(i-ml)b1-mb2 

xf (k+1) 

L  0 

p.A22+uA12 

xf(k) 

LB]_  +B2 

where  L  and  M  are  equilibrium  solutions  of 

^+1  =  ^LkA12Lk+  A22Lk"A21^  A11  Lo  *  "A21A11 
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\+l  =  All^A12LMk+MkA22 +MkLA12'A12^ 

A12  =  '^125  h.2  =  |iA22  Mo  =  ~A11A12 

and  A 

A  =  A,,  -  A,,L  has  N  slow  eigenvalues 
s  11  12  s 

Af  *  ^2  +  ^12  ^aS  Nf  ^ast  eigenvalues. 

For  p.  sufficiently  small,  a  designer  may  simply  ignore  ,  which 
is  equivalent  to  assuming 


l\i(Af)l  =  0  i  =  1 ,2 , . . .  ,Nf 

or,  all  our  fast  modes  are  assumed  to  have  a  deadbeat  response.  When 
this  perturbation  is  made  on  (4.11) 


x1(k+l) 
x2 (k+1) 


XL  (k) 

x2(k) 

(4.13) 


The  designer,  however,  uses  his  0(p.)  slow  model  and  chooses  some 
feedback  gain  G  to  place  the  Ns  slow  eigenvalues  of  the  reduced  order 
system 

x^(k+l)  =  A-qX^  (k)  +  B^u(k)  (4.14) 

where  u(k)  =  Gx^  (k)  providing  the  pair  (A^,B^)  is  controllable.  It  is 
interesting  to  note  that  (4.14)  will  result  if  p.  =  0  in  (3.32).  This  further 
verifies  that  to  0 (jj. ) ,  x^  (k)  will  determine  the  slow  dynamics  of  the  actual 
system. 

With  feedback  included,  the  closed  loop  system  is  now 


x ^  (k+1 )  (A11+B1G)  '  0  x^^) 

*2  (k+1)  •  (A21  +B2G)  |  0  x2(k) 


(4.15) 


which  preserves  the  deadbeat  approximation  of  the  fast  modes  in  accordance 


with  our  model . 
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How  accurate  is  this  reduced  order  design? 
i.  How  close  to  your  desired  "slow"  eigenvalues  are  your  actual 
eigenvalues  ? 

ii.  How  much  have  your  fast  eigenvalues  been  proved  as  a  result  of 
this  feedback? 

Lemma  8 :  Given  system  (4.5),  if  (A^,B^)  is  a  controllable  pair,  then  by 
choosing  Ng  slow  locations  for  eigenvalues  using  (4.11) 


'Vslow<'new)  =  Xslow('desired)  +  0(y'->  •  (4.13a) 

Xfast(new)  =*fast(0ld)  +0(^-  (4.13b) 

Proof :  Both  (4.13a)  and  (4.13b)  are  are  obvious  consequences  of  our  original 
transformation.  Our  original  system  (4.11)  with  this  reduced  feedback  looks 
like 


\  (fc+1) 

j 

au  +  big 

1 

■p 

1 

'  I-1 

x1(k) 

_x2(k+l) 

A21  +B2G 

1  ^^22j 

x2  (k) 

The  slow  eigenvalues  of  this  system  will  be  the  eigenvalues  of 
(An  +B1G)  +  H(A21  +B2G)(A12)(Au  +B1G)_1  +0(ii 2) 
verifying  (4.13a).  The  fast  eigenvalues  are  the  eigenvalues  of 

t*(A22-(A21  +B2G)(Au  +B1G)'1(A12))  +  00*2)  (4.15) 

as  compared  with  the  old  eigenvalues  of 


P- (A22 “A.2  i^i  iAf2 )  +  0  (p.  ) 


which  trivially  verifies  (4.13b),  but  more  importantly,  establishes  a  bound 


on  the  magnitude  of  the  fast  closed  loop  eigenvalues. 
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Example :  To  verify  these  accuracy  results,  let  us  try  varying  p  in  the 
following  system  after  we  generate  a  control  for  the  low  order  model. 
Given 


xL(k+l) 

_1  1  ! 

A 

xx  (k) 

1 

*2  (k+1 ) 

1  ^ 

0  1  , 

lo 

J 

x£  (k) 

1 

x3(k+l) 

1  * 

Ai 

1 

x3  (k) 

+ 

1 

_x4  (k+1) 

[.2  0  , 

V° 

1/ 

_x4  (k>_ 

1 

let  us  first  find  a  feedback  G  such  that 

u (k)  =  Gx1(k) 

and  the  eigenvalues  of  (A^  +15^)  =  .  96  +  .  lj  using  LINSYS  [16]  to  do  this 
we  get 

G  =  [ -.0116,  -.0684] . 

Now,  applying  this  feedback  to  our  original  system  we  will  show  the  actual 
eigenvalues  for  varying  values  of  u . 


4 

Slow  Eigenvalues 

Fast  Eigenvalues 

0 

.96  +  .lj 

0,  0 

.01 

.96568  + .09582 j 

.009699,  -.001055 

.05 

.9874  + .0741j 

.04321,  -.04321 

.1 

1.0206,  1.0057 

.08617,  -.08479 

.4 

1.37,  .9248 

.3378,  -.27269 

If  the  0 (jj, )  design  is  unsatisfactory,  the  logical  approach  might  be  to  take 
another  term  in  the  series  of  (3.32)  and  thus  use  3.35)  given  here  again 


for  convenience. 
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xL(k+l) 

_ 

(AU  +UA12A21A11)  1 

^A12<-B2'A2 1A11B1^ 

x1(k) 

+ 

V 

z  (k+1) 

L  0  ! 

0 

x2  (k) 

I 

u(k).(4.17) 


By  using  this  model,  we  should  be  able  to  model  the  x^(k)  states 

2 

to  0(jj.  )  accuracy.  Thus,  the  question  arises;  to  what  accuracy  can  this 
model  be  used  to  place  the  slow  eigenvalues? 

2 

If  we  truncate  our  transformed  slow  variable  xg  (k+1)  to  0(p,  ) 
accuracy,  we  get 

Xg(k+1)  =  ^An-HJ.A12A21A^J]xs(k)  +  [B1-HiA"jA12(B2-A21A^B1)]u(k).  (4.18) 

We  immediately  notice  the  absence  of  the  delayed  control  term.  To  see  a 
reason  for  this,  let  us  look  at  the  components  of  (k)  for  structure  "V" 
using  the  "F"  transformation 

xs (k)  =  (I-ML)x1(k)  -  Mx2(k) . 

We  have  already  shown  that  ||m||  =0(p,).  Therefore 

xs(k)  =  x1(k)  +  0  Cj.  ) 
xs  (k)  #  x1(k)  +  0(p,2) . 

2 

Thus,  while  model  (4.17)  accurately  defines  the  x^(k)  states  to  0(p,  ),  we 
see  that  past  0(p.)  xs(k)^x^(k)  in  the  nonhomogeneous  case.  So,  for  eigen¬ 
value  placement  of  slow  modes,  it  is  necessary  to  consider  (4.18).  So,  to 
place  our  slow  eigenvalues  to  0(j*^)  accuracy  consider  the  feedback 

usOO  -  Gsxs(k) 


when  applied  to 


*s(k+l)  =  [A11-hJ.A12A21A1^]xs  (k)  +  [  (I-uA^  [a£2A2jA££B^  +  p.A^A^2B2]  u  (k) 


(4.20) 


and  design  G.  to  place  N,  slow  eigenvalues  of 


(A  +  B  G  )  . 

S  S  s  ' 


Then  we  may  form  our  control 


us(k) 


G  x  (k) 
s  s 


=  Gg[  (T-p. A^ 2A12A21A1 1 )x£  (k)  +tiA[[A1^x2  (k)] 

=  Gg  (I_^A^^A^2A2^A^^)x^  (k)  +^GsA^ApX1  (k) 

9 

and  place  our  slow  eigenvalues  to  within  0(p.~)  of  their  desired  values  while 
changing  our  fast  eigenvalues  by  0(u“).  Now,  what  if  the  discrete-time 
system  does  not  fit  structure  "V, "  but  does  fit  the  general  model?  When 
this  is  the  case,  we  will  show  how  the  same  reduced  order  model  result. 
However,  the  perturbations  must  be  carried  out  in  the  transformed  state 
domain . 


xg (k+1 ) 
x£(k+l) 


If  the  transformed  system  (2.43)  is  truncated  after  one  iteration 


AU-^(1‘j)Ai2MjA2iAiJ  1  0 

n 

*s  (k) 

~  -  1*  -1  7l-i)- 

0  1  p,A22“p.  A21A11^  Aii 

Xf(k) 

+ 

(I"Ai^(1’j)Ai2^ji2iAn)Bi+An^(1’j)Ai2B2 

U(k) 

lAilBl  +B2 

and  since  p.  is  a  positive  scalar, 
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us(k)  =  Gsxs(k) 

=  Gst(I-t*A^jA12A21A"J)x1(k)  +  u(1_j)A^A12x2(k)]  .  (4.25) 

So  as  4"*0,  x^(k)^xs(k)  in  general.  Thus,  your  actual  feedback  will  be 
of  the  form 

u(k)  =  GsXl(k)  +  GsM.(1"j)A11^12X2(k)- 

The  x0(k)  term  may  or  may  not  be  neglected  depending  on  the  value  of  j. 

2 

0  fa  )  and  higher  accuracy  reduced  order  models  are  resolved  in 


the  same  manner. 


Example  0(y  )  design:  Consider  the  system 


Xj,  (k+1) 

*2 (k+1) 
x3 (k+1) 
x4(k+l) 


.1  .5 

0  1 


xx  (k) 
x2  (k) 
x3(k) 
x  (k) 


+  u(k) 

1 


In  a  previous  example,  for  y.  =  .1,  .4,  our  0(y.)  design  proved  unstable 

2 

for  these  same  values  of  y. ,  let  us  try  an  0(y  )  design,  y.  =  .  1 


1.12 

.78 

1.299 

= 

x(k)  + 

0002 

.  9998__ 

[_1 . 001 

u(k)  . 


For  \ ' s  of  . 96  + . lj 


ug(k)  =  (-. 038122,  -.15013)xgCk) 

=  (-.033525,  -.15854)x1(k)  +  (-.0038122,  -0039242)x, (k) 

which  gives  us  closed  loop  eigenvalues  of 

.95534  +  . 10174 j ,  -.086308,  .085836 

opposed  to  open  loop  eigenvalues  of 

1.1238,  .98549,  -.09298,  .084659. 

For  y.  =  .4 


x  (k+1) 
s 


1.48 

.12 

x  (k)  + 

2.16 

.008 

.992_ 

s 

.1 . 04_ 

us  (k) 


(-.23814  -.036177)xs(k) 


*  (-.12545  - . 24412 )xx  (k)  +  (-.095256 


- .  087177)x7 (k) 


which  gives  closed-loop  eigenvalues  of 
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.97222  +  . 2019 j ,  -.3854,  +.32901 

opposed  to  open-loop  eigenvalues  of 

1.4,  .98689,  .3317,  -.2786. 


4.5.  High -Accuracy  Reduced  Order  Modeling  for  Eigenvalue  Placement 

In  the  previous  section,  we  have  shown  how  reduced  order  models 
can  be  used  to  place  Ng  "slow"  eigenvalues  to  0 (jj. )  and  0 <jj.  )  accuracy. 
However,  what  if  for  any  one  of  a  number  of  reasons,  you  would  want  to  be 
able  to  place  the  fast  eigenvalues  at  fast  stable  locations.  By 

approaching  the  problem  from  our  limiting  cases  "V"  and  "H"  first,  we  hope 
to  point  out  the  advantages  and  disadvantages  that  each  transformation 
possesses  before  applying  both  transformations  to  the  control  design  of  our 
general  model.  Also,  what  if  our  control  algorithm  involved  an  adaptive/ 
identification  procedure  such  that  gains  had  to  be  updated  in  some  multi¬ 
processing  environment.  It  would  then  be  desirable  to  have  an  algorithm  of 
low  order  computationally,  but  very  accurate  that  could  be  easily  updated 
if  necessary. 

In  this  section  we  will  propose  two  high-accuracy  reduced-order 
modeling  algorithms  for  the  purpose  of  eigenvalue  placement.  Again,  the 
algorithms  will  be  dual  to  one  another.  Also,  various  structural  properties 
of  the  system  will  enable  us  to  improve  the  algorithms  efficiency. 

We  will  first  consider  an  algorithm  for  structure  "V."  A  dual 
procedure  for  "H"  will  be  outlined  in  much  less  detail.  Finally,  an 
algorithm  for  our  general  model  will  then  be  given  in  our  next  section. 
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Part  I:  Structure  "V1 


Given 


^11  I  ^12  x^(k) 

- =  -  ~  r  -  - +  —  u  (k ) 

x2(k+l)  A,,,  ,  p,A„„  x,  (k)  B„ 


M,a22J  [x2  (k)j  B2 


where 


i.  (k)  €  R  ,  x_2  (k)  €  R 

ii.  A^ j ,  Bj ,  i,j=l,2  are  real  constant  coefficient  matrices 

iii.  jj,  is  a  small  positive  parameter  such  that 


0<p.  < 


Iterate  the  following 


3(||aJJ||(||aoIH|lo||||a12H)  A°  kl2  L°Al2 


Lr+1  ^'k22\  +M'Lk^l2Lk'A21^All  Lo_A21All' 


(4.26) 


The  accuracy  you  desire  will  depend  on  the  number  of  iterations.  Letting 
"L"  be  the  equilibrium  value,  the  kth  iteration  will  give  the  following 


accuracy 


k+1 

L  =  1^  +  0(|/  V 


Now,  place  fast  modes  using  the  fast  subsystem 


xf(k+l)  =  M.[A22 +LA12lxf  (k)  +  [LB1+B2lu(k) 


(4.27) 


uf(k)  =  Gfxf(k) 


=  Gf[Lx1  (k)  +x2  (k)] 

=  GjLXj^  (k)  +  Gfx2  (k) 


(4.28) 


It  is  important  now  to  look  at  this  partial  feedback  in  both  our  transformed 
system  and  original  system.  In  the  transformed  system, 
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x  2  (k+1 ) 

*uV‘^  i 

xx  (k) 

+ 

V 

x2 (k+1) 

LA2i+U„B*G*fL  1  >*M«22+BIGPj 

x2  (k) 

.  B2  . 

uOO 


where 


Now,  if 


i*  = 
*12 


°<^M< 


_ 

*11  _ 

1  P‘MF12 

xx(k) 

+ 

V 

u(k) 

_F21 

u  F 

1  HyT22_ 

x2  (k) 

B2 

L  J 

A12^ 

-A. 

BlGf  = 

B. 

lGf 

B2 

— 

B2Gf 

u  * 

Gf 

(A. 31) 


3||f"Jii  (Hfoii-hil*iii!f12ii  > 


(4.32) 


Then  the  following  iteration  will  converge 


-1 


K+i  =  [^V^m^2^-f21]fu 

-1 


(4.33) 


Lo  " 


'F21F11' 


We  now  state  a  lemma  which  relates  L  to  L*  and  could  save  computations. 
Lemma  9 :  If  >Ng  and  if_  the  matrix  (4A.^2+^^G^)  is  of  full  rank,  the  L 
matrix  is  invariant  to  the  partial  feedback 

uf(k)  =  GfLx1(k)  +  G^  (k) 

and  therefore  L  =  L*. 

Proof :  As  we  have  already  shown,  the  partial  feedback  will  not  affect  the 
structure  of  the  homogeneous  s low  submatrix.  Therefore 
(F^  l"M,j^f,12^'  )  -  (<li  i  ^L) 


ll  12 

*  -t-n*r*\T*  = 


(An+P^B^GfL)  "  4m(A^2+B^Gj)L  A^  “  ^A]_2L 

A11  +tJ,MBlGfL  "  M'MAT2L'  '  ^MBlGfL  =  All_lJ,A12L' 
However,  since  ^{^-*2  =  ^A12  ’  *  BlGf 


■■■■■ 
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B,  G-L  -  P-A,  0L*  -  B,  G,L*  =  -HA.„L 


B1Gf(L-L*)  +  lJ*A12  (L-L*)  =  0 


OAi2+BiGf)(L-L*)  =  0 


Matrix  D  is  NgxN^.  Therefore  to  guarantee  (L-L*)  £  (D) ,  D  must  be  such 

that  N  >N  and  is  of  full  rank.  With  these  conditions  imposed 
f  s 

L  *  L*.  (7j  £  Null  Space) 

Now,  if  it  is  also  desired  to  place  the  slow  eigenvalues  to  a  desired 
accuracy,  you  would  now  iterate 


"k+l  ■'iMFU<f12“k+V22+Vi'l2-h2> 


Mo  *  -Vuhr 


Then  select  G  to  place  N  slow  modes  of  (A+B  G  )  where 
s  s  sss 


xs(k+l)  =  [An-p-A12L]xg(k)  +  [  (I-ML)B1-MB2lu(k). 


(4.34) 


Then, 


us (k)  =  Gsxs (k) 


Gg [  (I-ML)x1 (k)  -  Mx2  (k)] 
Gg  (I-ML)x1  (k)  -  GgMx2(k). 


(4.35) 


Our  composite  control  is  then 


uc  (k)  =  uf (k)  +  Ug  (k) 


=  GfLx1(k)  +  G^2  (k)  +  Gg(I-ML)x1(k)  -  GgMx,  (k) 
=  [GfL +Gf(I-ML)]xL(k)  +  [Gf-GgM]x2(k) 


(4.36) 


and  it  places  N  slow  \'s  of  (A  +B  G  )  and  Nf  fast  \'s  of  (Af+B,G,) . 

s  sss  I!  r  r  t 

The  design  has  been  carried  out  without  ever  handling  a  high  order  matrix. 
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Part  II:  Structure  "H1 


Given 


Xl(k+1)  Au  Au  Xl(k)  B 

- - - - - - +  -  -  u(k) 


x2(k+l)|  rA2i  |  ^A22  [_x2^  |_b2 


(4.37) 


where  i.  X^(k)€R  s,  j<2(k)€R 

ii.  A...  B.j  j=I,2  are  real  constant  coefficient  matrices 
ij  J 

iii.  u  is  a  small  positive  parameter  such  that 


0  <  M.  < 


sIIa-^i  (IIaJMIlJIIIaJ) 


Ao  =  A22_A21Lo ’ 


Then  we  can  iterate 


Hl-HI  A11(A12  +tJ‘LkA22'M'LkA21Lk)  Lo  A11A12 


(4.38) 


to  a  desired  accuracy. 


We  then  form  the  slow  subsystem 


xs  (k+1)  =  Asxs(k)  +  Bsus(k) 


where 


As=AU+^IA2i 

B  =  B,  +  LB- . 
s  1  2 

We  then  design  a  slow  feedback  Gg  to  place  the  Ns  slow  eigenvalues  of 

(A  +B  G  )  at  N  slow  locations, 
s  s  s '  s 

Our  partial  feedback  transformed  system  now  looks  like 


xs  (k+1 )  (a1I_Hj‘LA21+BsGs^ 


xf  (k+1) 


xs(k) 


V-  (A22+A21L)J  LXf  (k)  I  |-MB1+(I-ML)B2 


B1+LB1 


(4.39) 


If  this  system  were  reblockdiagonalized  through  an  successive  approxima¬ 
tion  iterative  procedure  (i.e.  x*(k)  = x^ (k) -K*xs (k) ) ,  our  designed  slow 
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eigenvalues  and  original  fast  eigenvalues  would  be  preserved.  The  only 
limitations  on  the  accuracy  would  be  how  many  iterations  of  (4.26)  we  used 
in  establishing  our  reduced  order  subsystem. 

When  implemented  on  our  original  model 


Us(k)  .  Gs,s(k) 

=  Gs[x1(k)  +Lx2(k)] 

=  G^Ck)  +  GsI^2(k)  (4.40) 

:x(k+l) 

2  (k+1 ) 

Depending  upon  the  shift  of  the  annulus  of  slow  eigenvalues  in  the  z-plane 
due  to  the  slow  feedback. 


All+BlGs  l  A12+B1GsL  Xl<k>  B1 


rA21+B2Gs  I  "A22+B2GsLJ  LX2^k>J  B2 


uf  (k). 


(4.41) 


!Im22+b2gsl||  . 

Ki+V.!l  '*  l”  sen'ral- 

Thus,  our  new  system  will  have  the  form 


At  this  stage,  the  slow  eigenvalues  have  been  placed  at  Ng  slow  locations 
to  a  desired  accuracy  using  a  reduced  order  slow  subsystem  (N  ) .  Plus,  the 
fast  eigenvalues  have  been  untouched. 

If  you  also  want  to  move  the  fast  modes  without  effecting  the 
slow  modes,  we  can  prove  a  dual  invariance  property  of  L  if  the  system  is  of 
the  proper  form.  To  block  diagonalize  system  (4.42)  two  matrix  iterations 
must  converge  and  thus  our  p  bounds  should  be  retested.  These  matrix 


iterations  are 
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Hc+l  =  FU^F12  +HcF22'LkF21Lk1 


(4.43) 


\+l  =  tF21+F22Mk+F21LX+MkL"F21]Fil- 


(4.44) 


Lemma  10:  If  N.  ^N-  and  A01  +B0G_  is  of  full  rank,  then  L  is  invariant  to 
'  SC  -L  _  s 

the  partial  feedback 

us(k)  =  GsXl(k)  +  GsLx2(k). 

Proof :  In  this  case,  the  fast  subsystem  does  not  change,  therefore 


F22  "  F21L*  ~  A22  “  ‘Si1 


A22  +  B2GsL  -  (A21  +B2Gs)L  =  ^2  ‘  A21L 
B2Gs[L-L*]  +  A21[L-L*]  =  o 

,[A21+B2Gs],[L~L*]  =  °* 

d(n^ns) 

Therefore,  if  >N.,  and  D  is  of  full  rank,  such  that 

[  L-L*] 

y- 

then  L  =  L  •  If  this  condition  holds,  we  would  then  iterate 


\+l  =  P‘^F21  +  F22Mk  +  F21LMk+MkLF21]Fir 


(4.45) 


Then,  using  either  (4.37)  and  (4.45)  or  (4.43)  and  (4.44)  we  would  place  N, 
fast  eigenvalues  using  the  fast  subsystem 


where 


xf(k+l)  =  Afxg(k)  +  3fuf(k) 
Af  -  uA22  -4A21Lk 


!f  *  -MB,  +  (I  +ML)Bn 


we  Chen  find  a  fast  control 


uf<k)  a  Gfxf(k) 


35 


Gf[  -Mx1(k)  +  (I  +ML)x2(k)] 
-GfMx1(k)  +  Gf  (I  +ML)x2(k) 


(4.46) 


x^(k+l) 

a 

x2  (k+1) 

A1 1  +3LG5'BlGrM’  A12  +B1GgL  +B1Gf  (I+ML) 
A21)+B2Gs'B2GfM’  ^2  AGsL*+B2Gf(I+ML) 


x2  (k) 
x,  (k)j 


(4.47) 


In  the  next  section,  we  will  give  an  example  based  on  model  "H."  In  the 
following  section  we  will  propose  an  algorithm  for  the  general  model. 

4.6.  Example  —  Structure  "H" 

Given 


x  (k+1 ) 


1  .6  |  .4  .8~ 

1 

-.2  1.3  1  0  .7 

1 

x  (k)  + 

— 

0  .1  1  .05  .7 

1 

-.1  .1  J  0  .1 

A 

u(k) 


The  open  loop  eigenvalues  are  at 


Therefore 


X.  ,  =  1.16  +  .37 j 

x3  =  .1 

X,  =  .0235. 


Nf  »  2 


The  desired  closed-loop  eigenvalues  are  at 


XL  =  .98 
X2  =  .96 

X3  =  *2 


*4  =  .15 


Iterating  the  L.^  matrix  to  6  iterations,  we  get 


H 
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.38416 

.45947 

.077  | 

. 65 1707  _ 

and  a  slow  subsystem  of 


1.045947  1 

.5156371 

1.8436 

= 

-  -  -  -1  -  —  - 

x  (k)  + 

— 

L-.  134829  i 

1.227129] 

s 

.1.7287 

us(k). 


Using  LINSYS  [16]  to  place  our  slow  modes,  we  found 


G  =  [ .077118  :  -.30607] 
s 

ug(k)  =  [. 077118,  -.30607,  .0060567,  -. 16403455] x (k) 


and  our  partially  closed-loop  system  looks  like 


x  (k+1) 


with  eigenvalues 


1.0771 

.29393 

.40606 

.63597 

-.12288 

.99393 

.0060567 

.53597 

.077118 

-.20607 

.056057 

-.064035 

-.022882 

-.20607 

.0060567 

-.064035 

X1  = 

.98005 

X2  = 

.95994 

U) 

II 

.099484 

II 

.023592. 

We  now  reiterate  both  the  and  matrix  recursions  to  show  the  invariance 
since  (Aj.^+I^Gs)  is  obviously  nonsingular. 


property  of 


L*  = 


M  = 


.38417 
_.  077 
f  .0605 


.45948 

.651713 

-.19712 


-.031422  -.138776 
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Our  fast  subsystem  is  defined  by 


xf  (k+1) 


.069815 

-.018603 


-.162899 

-.208847 


xf(k)  + 


1.4064 

1.5752 


uf(k). 


Again,  placing  our  desired  fast  poles 
Gf  =  [.64234,  -.42944] 

uf(k)  =  [-.048113,  .059422,  .62834,  -  .41281]  x(k) 
which  gives  closed  loop  eigenvalues  of 

=  .97999,  X  =  .96001,  X3  =  .19972,  X4  =  .15028. 


4.7.  Independent  Slow-Fast  Modeling  for  Reduced -Order  Eigenvalue  Placement 
In  the  previous  two  sections,  we  outlined  high-accuracy  reduced- 
order  modeling  algorithms  for  the  purpose  of  eigenvalue  placement  in  systems 
of  structures  "H"  and  "V."  While  we  were  able  to  prove  some  additional 
properties  of  certain  system  matrices  which  could  save  computations,  both 
algorithms  suffered  a  common  problem.  Given  a  system  fitting  model  "V," 
we  found  that,  in  general,  we  had  to  run  two  iterative  matrix  recursion 
equations  to  place  the  Ng  slow  modes  opposed  to  one  matrix  recursion  to 
place  the  fast  modes.  The  dual  problem  exists  for  model  "H."  This 
could  be  a  severe  time  setback  in  a  multiprocessing  control  environment. 

What  we  propose  now  is  to  apply  Section  2.5  of  Chapter  2  and  develop  an 
algorithm  similar  to  the  previous  section  such  that  it  can  be  applied  to  any 
system  of  our  general  form  (0^  j  i  1).  Also,  the  dependent  Q^(M^)  matrix 
recursions  will  never  be  needed. 
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3  A11  (  A22-A21A11A12  +  ^^11 


Likewise  for  transformation  "S"  to  converge 


3  ^1^  A22_A21A11A12  +  ^1  A11A12  ' 

Thus  the  logical  conclusion  since  each  bound  on  (4.49),  (4.50)  was  true 
for  any  j  6 [0,1]  .  Now,  for  any  satisfying  Lemma  ll,  we  propose  the 
following  straightforward  procedure.  Given 


xx(k+l) 


An  I  14  a12 


x2(k+l)  J  ^^2  J 


x,  00 


N  N,  M 

where  i.  x^(k)€R  ,  XjOO^R  ,  u(k)6R 

ii.  (A,B)  is  a  strongly  controllable  pair  in  ^ 

iii.  A^,  Bj  are  real  constant  coefficient  matrices 

iv.  0  £  j  £  1 


v.  satisfies  Lemma  11 


vi.  A, 


(1-1)-  -  1- 
U-'  A12’  ^22  a  ^22’  ^1  =  ^  A21‘ 


Then  define  a  slow  subsystem 
As  "  A11  +  PkA21 


Pk+1  "  ^1^12  +PkA22'PkA2lPk)  Po  *  A11A12 


with  control  matrix 


B,  "  B1  +  PkB2 


such  that 


x8(k+l)  -  Agxg(k)  +  Bgug(k) 


(4.52) 


and  define  a  fast  subsystem. 


Af  "  ^2  +  A12Pk 
a+1  "  (A22Pk  +  PkA12Pk'A21)All 


with  control  matrix 


Bf  -  PkBl  +  B2 


such  that 


Po  “  ‘*21*11 


xf(k+l)  -  AfXf(k)  +  Bfuf(k), 


The  strong  controllability  condition  in  p.  on  our  original  system 
guarantees  controllability  of  the  pairs 

(Ag,Bg),  (Af,Bf). 

The  designer  now  has  the  option  as  to  which  subsystem  he  wants  to  work 
with.  Let  us  say,  for  example,  that  the  designer  first  wishes  to  optimize 
or  place  the  slow  modes.  Then,  we  find  a  slow  feedback 

usOO  -  Ggxg  (k) 

for  (4.52),  where  G  is  the  product  of  an  N  order  pole -placement  procedure 
s  s 

or  a  slow  quadratic  performance  criteria 


WV  ■  4k£0C«;oOQ,*,W+u'(k)RuJ(k)> 
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where  (A,D  )  is  an  observable  pair,  D 'D  ■  Q 
s  s  s  s 

R  >  0  Q  >  0 
s 

Ug(k)  -  Ggxs(k) 

=  -(R+B^K  Bs)'1(B^KAg)xs(k) 

where  K  is  the  solution  to 

K  +  %  ■  <AsKV  «  +  «aKss)1  (sjK*,)  • 

Either  method  of  finding  G  is  acceptable  as  long  as  u  (k)  *  G  x  (k)  places 

s  s  s  s 

the  N  slow  modes  at  N  slow  locations  (i.e.  the  slow  modes  remain  in  some 
s  s 

stable  annulus  sufficiently  spread  from  the  fast  modes  which  are  within  some 
radius  R(p.)  of  the  origin  in  the  z-plane).  This  requirement  is  necessary 
if  we  are  to  next  carry  out  a  reduced  order  fast  design. 

Once  Gs  has  been  obtained,  to  implement  this  control  on  the  original 

system, 

us(k)  -  Gsxs(k) 

-  Gs[*1Ck)  +Psx2(k)] 

-  tGs  •  GgPg] x (k) .  (4.53) 

This  simple  form  of  the  solution  enables  our  control  to  span  only 
the  slow  controllable  subspace,  and  thus  keep  all  fast  modes  invariant  to 
0  Ox  ).  This  is  seen  much  in  the  same  way  as  in  the  previous  chapter.  Our 
"S"  transformation,  when  applied  to  our  general  model,  results  in 


~xg  (k+1) 

A  !  0 

s  1 

*s  00 

+ 

B 

s 

xf  (k+1) 

1  i 

°i 

1 _ 

_Xf(k) 

■  B"f 

u(k) 
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where 


A11  +  PkA2l 


B1  +  PkB2 


Af  “  *22  ~*21Pk 

Bf  "  +  (X'<Pk>B2 


with  the  feedback  u  (k)  ■  G  x  (k)  applied 

s  s  s 


<V*.G.>  . 


xs  (k+1)' 

_xf(k+!)_ 

BfGs 


xs  (k) 

+ 

V 

xf  (k) 

Bf. 

u(k) 


(4.54) 


We  could  block  diagonalize  this  system  by  letting  (u(k)*0) 


xf (k)  -  xf (k)  -  Txg(k) 


x*  (k+1 )  ■  x^(k+l)  -  Txg  (k+l ) 


(BfGg-T(Ag+BgGg)  +  AfT]xg(k)  +  A^Jfk) 


we  want  T  such  that 


BfGg  "  T<AS  +  BsGs)  +AfT  *  °* 


(4.55) 


However,  since  (Ag+BgGg)  and  A^  have  no  eigenvalues  in  common,  the 

existence  of  a  solution  to  this  Lyapunov  type  equation  is  guaranteed  by  [14]. 

Thus,  the  eigenvalues  of  (4.54)  will  be  determined  by  the  eigenvalues  of 

(A  +B  G  )  and  A-  preserving  the  accuracy  of  our  designed  slow  eigenvalues 
s  s  s  c 

tc 

and  the  original  fast  eigenvalues  to  OGi  ). 

If  the  designer  now  wants  to  generate  a  reduced  order  fast  sub¬ 
system  for  design  purposes,  he  must  first  update  the  system  matrices 


*11  *  A11  +  BlGs’ 

A* 

a12 

*21  *  *21  +  B2Gs  ’ 

CM 

+  BlGsP 

+  b2g3ps 


(4.56) 
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mm 


Since  ||p£||  =  OO*1"^),  if  ||B1Gg||  *  0(1)  and  IIb^JI  =  0(p^),  or  less,  our 
bounds  on  will  be  unchanged.  The  bound  on  ||b,G  ||  is  easily  satisfied  if 
we  remember  the  set  of  eigenvalues  we  are  considering  are  within  some  radius 
1+e  of  the  origin  of  the  z-plane.  The  bound  on  lU-G  II  is  somewhat  limiting 
and  suggests  the  general  form  B2=M*^B2.  However,  properties  like  this  come 
down  to  a  case  by  case  study.  In  general,  after  the  updates  (4.56),  define 
the  fast  subsystem 


a£2  +  A*„P,f 


“f  "22  12  k 

PL+1  *  (A22PL  +  PLA12PL_A21)A11  Po 
with  control  matrix 


"*21*11 


-1 


Bf  =  PlBi  +  B2 


such  that 


xf(k+l)  *  Afxf(k)  +  Bfuf (k) , 


(4.57) 


Find  a  fast  feedback  matrix  through  pole  placement  or  a  fast  order 
optimal  linear  regulator  solution.  Either  way,  we  have  fast  modes  at 
new  fast  locations.  To  implement  on  actual  system, 

uf00  “  Gfxf(k) 

■  Gft  P^xt(k) +x2(k)] 

“  [GfpJ  :  Gf]  x  (k)  . 

Again,  this  simple  form  of  the  solution  enables  our’ control* to  span  only  — 
the  fast  controllable  subspace,  and  thus  keep  all  slow  modes  invariant  to 
0&tL). 

Combining  this  with  our  previously  defined  slow  control  we  get 


u  (k)  *  ug(k)  +  uf(k) 

-  (Gg+Gfpf  ;  Gf  +  GgP®]x(k) 


(4.58) 
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which  places  N  slow  poles  at 
s 

Ua3+b9gs)  +  oo*k) 

and  fast  poles  at 

\(Af +  BfGf)  +  0(p.L). 

The  order  of  fast  and  slow  designs  is  not  important  as  long  as  the  dual 
procedure  is  carried  out.  Again  note  that  we  could  have  specified  our 
control  matrices  to  be  of  some  special  form  (i.e.  ^  ^B^)  and 

perhaps  simplify  the  computations.  Instead,  by  imposing  the  less  conserva¬ 
tive  form,  HbJI  ■  0(1)  or  less,  we  encompass  a  greater  class  of  systems 
possessing  open- loop  two -time -scale  properties. 

In  the  next  section,  we  give  an  8th  order  power  system  example. 

A  flow  chart  Illustrating  this  design  procedure  is  given  on  the  next  page. 


with  open  loop  eigenvalues  of 


Check  Bounds  on  fj. 


Figure  4.1.  Flowchart  for  independent  slow/fast  subsystem  design. 


\1  -  -5.0066 


\2  “  -4.3542 
\3  -  -1.9939 
\4  -  -1.3884 


X5  -  -.086912  +  .33047 j 

-  -.086912  -  .33047 j 

-  -.37062 
X8  -  -.075716 


which  would  suggest  N  ■  Nf  *4. 

s  c 

This  system  is  then  run  through  a  continuous  singularly  perturbed 
modeling  software  package  in  an  attempt  to  put  this  system  into  continuous 
two -time -scale  form.  As  a  result,  the  above  system  is  transformed  by  the 
permutation  matrix 

P  “  [e1,e3>e6,e8,e5,e4,e7,e2] . 

This  system  is  to  be  controlled  in  a  sampled-datn  environment 
where  Is  .9  sec.  The  control  design  is  then  best  carried  out  in  the 
z-plane.  The  corresponding  discrete  model  is  given  as 


x  (k+1) 


.83527 

0  0 

0 

0 

0 

0 

0 

.096874 

.86071  0 

0 

0 

0 

0 

.0293 

-.00223 

-.005  .88241 

-.25361 

.040567 

-.0025737 

-.0251 

-.00111 

.0068261 

.014284  -.029647 

.92793 

0 

.00576 

.059 

.00227 

-.029713 

-  .0609  2.0283 

-2.303 

.08842 

-.021163 

- .22386 

-.008 

.047959 

.75863  0 

00 

0 

.1653 

0 

.022615 

-.012274 

-.02676  1.2091 

-1.4005 

.16097 

-.01326 

.15593 

-.00558 

_.  81558 

0  0 

0 

0 

0 

0 

.01111 

3.2946 

00 

.15247 

-.0025695 

.010316 

-.050766 

.056549 

-.01494 

2.4774 


00 

.2936 

-.038309 

2.7624 

00 

1.4729 

0 


uOO 


(4.59) 


r 
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with  system  eigenvalues 


.874  +  .  1692  j 

X  -  .2866 

.874  -  . 1692 j 

X6  «  .1653 

.8607 

\7  =  .02 

.8353 

Xg  -  .0111 

Thus,  our  two-time-scale  eigenvalue  separation  is  preserved. 

After  testing  (4.59)  with  some  design  software  [Appendix  B] ,  we 
found  (4.59)  to  be  of  general  form  with  =  .259904.  If  a  reduced  order 
slow  model  is  desired,  transformation  "S"  is  generated  using  2  iterations 


of  Pk  giving 


0 

0 

.048604 


0  0  0 

0  0  .0344832 

0011815  -.0075897  -.0003338 

.0028456 


.0118520  .0072417  .0877587 

3 


Since  0(p.  )  accuracy  will  be  sufficient  for  slow  modes.  We 
thus  generate  the  fourth  order  model 


xs  (k+1) 


.83527 
’  .125 
-.0039 
.008065 


0  0 

.8607  0 

.008567  .971817 

.01671  .1005 


0 

0 

■.3549061 

.777731 


*soo 


3.2946 

.2379 

-.00582 

.01586 


0 

0 

.416685 

.123691 


UsOO 


(4.60) 


with  eigenvalues 
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•  • 

4* 


I 


\l  =  .87477  +  . 16202 j 

\2  -  .87477  -  . 16202 j 

\3  -  .80071 

-  .83527. 

4 


The  designer  now  wishes  to  move  the  open  loop  slow  modes  to 
(.9  +  .05J,  .85,  .80).  Using  LINSYS  [16]  on  (4.60),  we  obtain  a  feedback 


gain 


.075043  -.60196  0  0 

-.073167  -.036092  -.46533  .76279 


and  when  implemented  on  our  actual  system  this  becomes 


ug(k)  -  Gsxs(k) 


Gg^Ck)  +P®x2(k)] 


[G  •  G  Pf]x 


[- 


s  •  s  1  - 
.075043  0.60196 

-.073167  -.036092 

0 


0  0  ; 

•.46533  .76279', 

0  0 


,-.013577  .0060737  .070473 

and  gives  closed  loop  eigenvalues  of 


\l  “  .9 


.05  j 
.9  +  .05 j 
.80976 
.84142 


-.020758 

.0010813 


8 


.28400 

.1653 

.019959 

.011109 


which  verifies  the  0(p.  )  accuracy..  Now,  if  the  designer  also  chooses  to 
move  the  fast  modes,  the  system  submatrices  are  updated  and  the  "F" 


transformation 

is  updated 

using  3 

iterations  of 

_f  .  . 

Pk  giving 

.143835 

.5253375 

-.988291 

.181422 

P3f  “ 

.066625  - 

.059635 

•.9768864 

.45438 

0 

-1.3080671 

0 

.4155951 

-.896520  - 

.378294 

0 

0 

and  a  fourth  order  fast  subsystem  model 


x  f  (k+1 ) 


0148718 

-.00261 

-.0146517 

-.00164 

0 

.1653 

0 

-.00864 

0933589 

-.001 

.28884 

.004763 

0 

0 

0 

.01111 

+ 


.50762 

.1271 

.2584564 

-.533954 


2.4653 

0 

1.0729305 

0 


uf(k) 


with  eigenvalues  at 

=  .01996 
\2  =  .28375 
X3  =  .16530 
\4  -  .01111. 


The  designer  wishes  to  place  these  eigenvalues  at  (.2,  .1,  .05,  .01). 
LINSYS  [16]  on  (4.61)  we  obtain  the  feedback 


Gf  * 

0 

.09359 

0 

-.20918 

t 

_-  .044114 

. 16405 

-.12578 

-.011122 

When  implemented 

on  actual 

system 

becomes 

(4.61) 


Using 
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uf(k)  =  G^jCk) 

=  Gf[  PfxL  (k) +x2  (k)] 

=  [GfPf  j  Gf]  x  (k) 

’.19377  -.01229  0  0 

_.  0070544  -.23638  .20813  -.060278 

|  0  .09359  0 

I 

• -.044114  .16405  -.12578 

Thus,  the  net  control 

u(k)  *  ug(k)  +  uf(k) 

results  in  closed  loop  eigenvalues 

*  .89999  +  ,050005j  = 

\2  -  .89999  -  . 050005 j  X&  = 

X3  *  .84268  X  = 

X.  *  .80758  \Q  = 

4  8 


I 

I 

I 

I 


i 

i 

i 

I 

-.20918  " 

x(k). 

-.011122 


20000 

099952 

056108 

01101. 
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5.  CONCLUSION 

la  this  thesis  the  concept  of  a  two-time-scale  linear  shift- 
invariant  discrete-time  system  is  introduced.  Properties  of  such  systems 
have  been  well  documented  for  the  continuous  case.  These  properties  include 
reduced  order  modeling,  separate  slow-fast  control  law  design,  boundry 
layer  phenomena,  etc.  It  was  the  purpose  of  this  thesis  to  develop  a  model 
for  discrete- time  systems  which  possessed  as  many  of  these  properties  as 
possible  without  restricting  the  model  to  isolated  special  cases. 

The  main  results  will  now  be  reviewed  very  briefly.  First, 
criteria  for  the  two-time-scale  property  of  a  system  using  matrix  norm 
conditions  are  proposed.  Then  a  model  satisfying  these  conditions  is  intro¬ 
duced.  If  a  two-time-scale  system  matrix  does  not  fit  the  model,  a  simple 
procedure  of  permutation  and  scaling  is  outlined  such  that  the  transformed 
system  matrix  will  fit  the  model.  A.  pair  of  dual  transformations  are  pro¬ 
posed  which  transform  our  general  model  into  slow-fast  block  diagonal  form. 
Conditions  are  given  for  the  existence  of  both  transformations.  The  difference 
between  a  degenerate  and  normal  system  is  briefly  discussed. 

The  recursive  property  of  discrete-time  systems  is  used  to  reveal 
the  two-time-scale  nature  of  our  model  by  iteratively  solving  for  slow  or 
"steady-state"  models.  We  thus  introduce  reduced  order  modeling  and  the 
concept  of  a  "boundry  layer"  for  our  class  of  systems. 

Separate  slow  and  fast  control  law  design  is  then  established  for 
the  purpose  of  eigenvalue  placement.  Controllability  conditions  are  dis¬ 
cussed  and  an  algorithm  is  proposed  for  reduced-order  control  law  design 
which  maintains  is  high  level  of  accuracy  when  implemented  on  the  high 


101 


order  system.  Also,  the  basic  differences  between  discrete  two-time-scale 
systems  and  continuous  singularly  perturbed  systems  are  discussed. 


H 


APPENDIX  A 


PROOF  OF  CONDITIONS  NEEDED  FOR  BLOCK  DIAGONALIZATION 


The  proof  given  in  [6]  is  applied  several  times  in  this  thesis. 
For  the  convenience  of  the  reader  it  is  restated  here  in  context. 

It  is  required  that  L  exists  such  that 


A22L  -  LA11  +  LAi2l  *  A21  *° 


L0”A22A21  Ao“A11’A12L0 


Bl*All'A12L  B2*A22  +  LA12 


seek  L  of  the  form  L*Lq  +  D.  Where  D  is  a  real  root  of 

DA0-(A22  +  LoAi2)I>-n*l2D  +  '-o*0a0-  « 

The  following  lemma  gives  a  sufficient  condition  for  the  existence  and 
uniqueness  of  a  real  root  D  and  establishes  a  bound  for  its  norm  j|  D  ||  . 

It  also  formulates  a  convergent  procedure  for  iterative  calculation  of  D. 
Lemma  1:  If  A22  is  nonsingular  and  if 


!Ia22II  <-j|Iao!I  +  IIai2!IIIloI 


then  a  unique  real  root  of  (A2)  exists  satisfying 


0 <  UD!|  <■ 


2llAo!!llLoi _ 


This  root  is  an  asymptotically  stable  equilibrium  of  the  difference  equation 


Dk+1  ■  *22  ‘Vo  +  DkA0  •  ■  Dk412Dk>  *f  <Dk> 
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and  its  domain  of  attraction  encompasses  the  set  of  matrices  defined  by 
(A4). 

Proof:  For  a  —  IWqII*  b  ”  !IA12^^L0^*  c  "  ^22'^  and 

IIdJI 


IAqII  l'L0 


(A6) 


we  obtain  from  (A5) 


dk+1<c[l+(a  +  b)dk  +  abd^<cj^  dk+  l) 


(A7) 


and  we  analyze  the  upper  bound  of  d^  defined  by 


dk+l  ’ ' 


a  +  b  -7  . 

~  dk+1 


(A8) 


Obviously  dk<dk  for  all  k.  When  c<  j(a  +  b)  ^  then  the  scalar  difference 
equation  (A8)  has  two  real  equilibrium  points  d'  and  d" >d'  and  d'  lies  in 
the  interval  [0,  2 (a  +  b)  .  For  all  d^^d'  in  this  interval  we  have 
ldk+^  -  d '  I  <  Id^  -  d' I ,  which  proves  that  is  bounded,  that  is 

0  <||dki|<2(a  +  b)"L  (A9) 

holds  for  all  k*l,2,...,  if  it  holds  for  k  ■  0.  Substituting  ^Dk*Dk-D 
into  (A5)  we  get 

50k+l  *A22l6Dk<A0  -  A12D)  -  (L0  +  D)A126Dk  -  6DkA126Dk]  (A10) 
which,  using  |i6Dk||  « vk<ll  D^ll  +  lloll  »  i®plies 


\+i  <cta  +  b  +  3||a12||||d!I+||a12|||| dJ|  ]  vfc. 


(A1 
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When  both  ||dJ1  and  ||d||  satisfy  (A9) ,  then  (All)  yields 


Vi<c(a+b+  !rb)vk<3e<*+l,)',k 


(A  12) 


and,  hence,  if  c<-|(a  +  b)  \  as  required  by  (A3),  then  f(D^)  in  (A5)  is  a 
contraction  mapping  and  D  is  its  fixed  point.  We  complete  the  proof  by 
noting  that  is  a  Lyapunov  function  and  <  0  f°r  D(c  satisfying 

(A9) . 

Using  (A5)  with  the  initial  condition  Dq = 0  we  can  calculate  D 
iteratively.  By  (A3)  and  (A12)  after  k  iterations  the  relative  error  is 


I1^11<[3||a2-21II<I1a0II+||AuIl011  )!k  (ai3) 

and  it  decreases  as  IIa^^II  and !l IIAq  11  decrease,  that  is  as  the  ill-conditioning 
of  the  system  increases.  When  we  need  M  to  be  a  real  root  of 


BjM  -MB2  +A12  -  0. 


(A  14) 


The  following  lemma  formulates  a  convergent  iterative  method  for  solving 
(A  14). 

Lemma  2:  Under  the  conditions  of  Lemma  1  the  solution  M  of  (A14)  is  the 
asymptotically  stable  equilibrium  of  the  linear  difference  equation 

Vl  ■  I  <AU  '  At2L>"k  -  MkLA121A22+A12A22'  <Al5> 

-11^  -Mil  we  obtain  from  (A14)  and  (A 15) 

mk+1<c[a  +  b  +  2|lA12ll  ||d||  ]v^  (A16) 


and,  by  virtue  of  (A4) 
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Thus, 


“k+l 


a  +  b  + 


4ab 
a  +  b 


V 


(A17) 


Vi<nk 


Cn^  if  c<  y(a  +  b)  1  which  is  satisfied  by  (A3) 
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APPENDIX  B 
COMPUTER  PACKAGES 

Through  the  use  of  the  IMSL  matrix  analysis  computer  library  and 
theA.G210  computer  graphics  library,  an  extensive  discrete  two- time-scale 
systems  software  packages  has  been  compiled.  This  software  accepts  an 
arbitrary  linear  shift-invariant  discrete-time  system  and  determines  whether 
or  not  it  possesses  a  two-time-scale  property,  what  order  reduction  is 
possible,  comparitive  time  responses,  reduced  order  eigenvalue  placement, 
etc.  Rather  than  list  the  source  program  here,  the  author  welcomes  any 
party  that  may  be  interested  to  contact  him. 
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